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Resumen 



La actual expansion acelerada del universe sigue siendo un gran misterio para la fisica. Varies modelos 
se han propuesto para su explicacion, entre ellos el de la energia oscura, que sin embargo posee proble- 
mas tanto de caractcr tcorico como experimental. Sin embargo, otras aproximaciones, como la de las 
teorfas de gravedad modificada son una interesante alternativa para resolver este problema. Motivados 
en este enfoque realizamos un estudio de los modelos cosmologicos en las teorfas f{R) que son extensio- 
nes naturales de la Relatividad General con funciones arbitrarias del escalar de Ricci. Un capitulo esta 
dedicado a la obtencion de las ecuaciones de campo en el formalismo metrico de las teorias f{R), inclu- 
yendo la discusion sobre los terminos de frontera. Aplicamos dcspues tales ecuaciones para escribir las 
ecuaciones dinamicas del universo, aceptando para esto la validez del modelo estandar de la Cosmologia, 
es dccir, usando como espacio-tiempo el de Robertson- Walker. La aproximacion de sistema dinamico 
nos permite encontrar criterios de viabilidad para las funciones f{R), cxigicndolcs que cumplan con una 
transicion de cpoca de dominio de materia, a una de expansion acelerada. Ccntramos gran atencion en 
el problema de las distancias cosmologicas en teorfas f{R), obteniendo expresiones generales a partir 
del estudio de la Ecuacion de Desvfo Geodesico. 
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Abstract 



The actual accelerated expansion of the universe continues being a mystery in physics. Some models 
had been proposed for this explanations, among them the dark energy, which however has problems of 
experimental character as well as theoretical. Other approximations, like modified gravity theories are 
an interesting alternative for this problem. Motivated in this approach we study cosmological models 
in f{R) theories which are natural extension of General Relativity with arbitrary functions of the 
Ricci scalar. One chapter has dedicated to obtain the modified field equations in the metric formalism 
of f{R) theories, including the discussion about boundary terms in the action. Later, we apply these 
equations in order to describe the dynamics of the universe, using for this as space-time, the Robertson- 
Walker metric. The Dynamical system approach for this equations allow us to find viability criteria for 
the functions f{R), demanding them to satisfy a transition between a mater dominated epoch to an 
accelerated dominating epoch. We focus our study in the problem of cosmological distances in f{R) 
theories, obtaining general expressions from the Geodesic Deviation Equation (GDE). 
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Notacion y Convenciones 



La convencion de signos usada para la metrica es (— , +, +, +), la velocidad dc la luz se toma en unidades 
geometrizadas c = 1. Los indices griegos a, P, . . . se utilizan para las coordenadas espacio-temporales y 
van dc a 3, los indices latinos a, 6, . . . se utilizan linicamente para las coordenadas espaciales y van de 
1 a 3. Letras en negrilla A, B, . . . son utilizados para denotar tensores. 



G Constante de Gravitacion Universal (6.67 x 10 m^/kg s^) 

K 8ttG 

Qap Componentes del tensor metrico 

g Determinante del tensor metrico 

Ra^-yS Tensor dc Ricmann 

Ral3 Tensor dc Ricci 

R Escalar de Ricci 

Simbolo de ChristofFel 

Gcep Tensor de Einstein 

da Derivada parcial con respecto a 

Va Derivada covariante 

A^B Derivada de A con respecto a B 

□ d'Alambertiano □ = V^V-y 

Derivada covariante a lo largo de la curva 

A(^a0) Simetrizacion 

A[Q,y3] Antisimetrizacion 

Derivada temporal 

Tajs Tensor de Energia-Momentum 

a{t) Factor de escala 

t Tiempo cosmico 

z Redshift 

H Parametro de Hubble 

k Curvatura espacial del universo 
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w Ecuacion de estado w = p/ p 

A Constante cosmologica 
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Capitulo A 

Prefacio 



"p^ esde el nacimiento de la Relatividad General, que se ha proclamado como la teoria de gravitacion 
con mas aceptacion actualmente, se ha buscado una solucion que pueda describir la evolucion del 
universe. El modelo estandar de la cosmologia esta basado en dos principios fundamentales que nos 
dicen que a gran escala el universe es estadisticamente homogeneo e isotropico. La solucion de las ecua- 
ciones de campo de Einstein para un universo que satisface estas condiciones nos conducen a la solucion 
de Friedmann-Lamaitre- Robertson- Walker (FLRW) [1],[2]. 

Einstein intento conciliar el modelo de un universo dinamico introduciendo constante cosmologica A 
con el fin de tener un universo estatico. Sin embargo, desde las observaciones realizadas por Hubble en 
1929 que mostraban que las galaxias se alejaban mas rapido en cuanto mas lejos estuvieran, mostrando 
asi un universo en expansion, los cosmologos y la comunidad cientifica adoptaban cada vez mas la idea 
de un universo dinamico. Einstein considero entonces la introduccion de la constante cosmologica como 
el error mas grande de su vida. 



Hoy en dia sabemos que nuestro universo ademas de expandirse, lo hace de forma acelerada, esto 
gracias al analisis de los parametros cosmologicos a partir del estudio de Supernovas de tipo la (Sn 
la) [3]. Sin embargo no se ha encontrado una explicacion satisfactoria para este hecho. Por otro lado 
las observaciones del proyecto Wilkinson Microwave Anisotropy Probe (WMAP), en las anisotropias 
de la Radiacion Cosmica de Fondo (CMB), nos predicen que el universo esta compuesto en un 74% de 
una misteriosa forma de energia, de la cual no conocemos su naturaleza, y un restante 26 % de materia 
compuesta en un 22 % por materia oscura y el restante 4 % en materia barionica. Este tipo de energia 
lo denominamos energia oscura. Se pueden proponer 3 modelos para explicar la expansion acelerada del 
universo: 1. Retomar la idea de la constante cosmologica A, 2. Energia oscura, 3. Gravedad modificada 
[4]-[8]. 
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Capi'tulo 1. Prefacio 



El modelo de la constante cosmologica provee una explicacion para la expansion del universe, aun- 
que posee dos problemas que han Uevado a la comunidad cientifica a buscar nuevas alternativas. A 
partir de los datos observacionales se puede estimar un valor para la constante cosmologica, el cual 
difiere con valores obtenidos en teorias cuanticas de campos en 120 ordenes de magnitud; este problema 
es conocido como problema de la constante cosmologica. Otro problema fundamental es el denominado 
problema de coincidencia, y es el por que actualmente la densidad de energia de la constante cosmolo- 
gica es comparable a la de la materia [9]. 

Los modelos de energia oscura vienen de suponer una ecuacion de estado de la forma w = p/p, en 
donde p es la presion y /? la densidad de energia. Para un valor de w < — 1/3 se obtiene una expansion 
acelerada del universo, como lo muestran las ecuaciones de Friedmann; el caso w = —1 corresponde a 
un fluido de densidad constante que puede asociarse con la constante cosmologica [6]. Los modelos en 
los cuales — l<w<— l/3se denominan de Quintaesencia [10]. 

La tercera posibilidad consiste en ver los terminos adicionales en las ecuaciones de campo de Einstein, 
que dan la expansion acelerada del universo, a modificaciones de la gravedad. Una de est as formulaciones 
son las denominadas teorias de gravedad modificada f{R) que surgen de una generalizacion de la accion 
de Einstein- Hilbert con terminos adicionales del escalar de curvatura R [11]- [7]. Existen tres versiones 
de la gravedad f{R)- 1- Formalismo metrico, 2. Formalismo de Palatini y 3. Formalismo metrico-afm; 
estos tres formalismos se diferencian basicamente en las dependencias de los campos de materia con el 
tensor metrico y las conexiones, y tambien en la existencia de una conexion afm en la variedad. 

Puesto que estas teorias se pueden ver como extensiones de la Relatividad General, deben producir 
los mismos resultados, por ejemplo tener limites Newtonianos correctos a nivel del sistema solar, poseer 
un problema de Cauchy bien definido y no sufrir de inestabilidades (en el sector de materia por ejemplo 
existe la llamada inestabilidad de Dolgow-Kawasaki [18], existente en algunos modelos propuestos [19]; 
por otro lado algunos modelos muestran que el universo se acelera por siempre o colapsa en un tiempo 
finito, las cuales se estudian en la denominada inestabilidad de De Sitter [20]). Existen varias formas 
de la gravedad f{R) que predicen de forma acertada la epoca de inflacion y la expansion acelerada del 
universo [21]. 

El objetivo principal del trabajo es obtener modelos cosmologicos en gravedad modificada f{R) que 
nos den cuenta de la expansion acelerada del universo, y ademas centrarnos en el problema de las 
distancias cosmologicas en gravedad f{R), que es el punto central de la tesis. Se estudiara ademas la 
teoria de perturbaciones cosmologicas (a primer orden) y el problema de Cauchy en teorias f{R)- Para 
esto supondremos de nuevo que el universo es homogeneo e isotropico, y estudiaremos las ecuaciones 
de Friedmann modificadas en este formalismo. 
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Capitulo 
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Elementos de Relatividad General 



'Most of the fundamental ideas of science are essentially 
simple, and may, as a rule, be expressed in a 
language comprehensible to everyone." 

Albert Einstein. 



a Relatividad General (RG) es la teon'a de gravedad mas aceptada actualmente. Propuesta por 



Albert Einstein en 1916, la Relatividad General explica los fenomenos gravitacionales como conse- 
cuencias de la curvatura en el espacio-tiempo generada por el contenido de materia presente. De esta 
manera se pudo realizar una transicion de la descripcion fenomenologica de la gravedad vista como una 
fuerza de atraccion (Newton) a una propiedad netamente geometrica del espacio-tiempo (Einstein). 
Aunque la RG cuenta con un fundamento observacional que la hace una de las teorias mas solidas en 
fisica, la gravedad sigue siendo aun una de las interacciones mas enigmaticas. Nuestro proposito en este 
capitulo es introducir muy rapidamcnte los postulados y ecuaciones de la RG, haciendo sin embargo 
gran enfasis en los principios variacionales de la RG, las cuales seran la base fundamental de la tesis. 

2.1. Principios de Relatividad General 

Como bien es sabido la gravedad fue primeramente descrita por Newton como una fuerza que actiia a 
distancia y que depende principalmente de las masas de los objetos. Esta descripcion permitio reali- 
zar grandes avances en la fisica, especialmente en la astronomia y condujo a valiosos descubrimientos. 
Aunque la descripcion Newtoniana de la gravedad permitia explicar fenomenos que tiasta el siglo XVII 
habian permanecido son resolver (en especial el movimiento de los astros en el sistema solar, sin contar 
la anomalia en la orbita de Mercurio, y la caida de los cuerpos), presentaba algunas fallas, en especial 
la imposibilidad de explicar el origen de esa fuerza de atraccion inherente a todo cuerpo con masa. 

La teon'a Newtoniana basa sus leyes en la existencia de los sistemas de referencia inerciales (siendo 
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2.2. Postulados y Formulacion Matematica de la Relatividad General 
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sin embargo la existencia de tales sistemas privilegiados un problema actual en fisica), tiene un gran 
problema cuando su fuerza de atraccion es considerada, pues presupone una interaccion a distancia 
con una velocidad de informacion infinita. Esto llevo a Einstein en 1905 a postular la Teoria de la 
Relatividad Especial (RE) basada en los siguientes postulados [22]: 

1. Las leyes de la fisica son independientes del sistema de referenda inercial. 

2. La velocidad de la luz c en el vacio es la misma para todos los observadores inerciales, indepen- 
diente de la direccion de propagacion. 

El caracter fundamental se basa entonces en la constancia de la velocidad de la luz, razon por la cual 
la teoria gravitacional de Newton deben'a ser reformulada. En esto trabajo Einstein entre los aiios 1905 
a 1916 en su Teoria de la Relatividad General (RG). Una de sus bases fundamentales es el denominado 
Principio de Equivalencia de Einstein que generaliza el Principio de Equivalencia Debil ^ y que puede 
escribirse como: 

Principio de Equivalencia de Einstein: No es posible determinar bajo ningun experimento si un objeto 
se encuentra en caida libre en un campo gravitacional uniforme 5 o si esta en un movimiento uniforme- 
mente acelerado con aceleracion a = g. 

En donde g es la aceleracion del campo gravitacional. La idea de Einstein fue entonces ver la tra- 
yectoria de partfculas en caida libre y ver a que efecto se debe su aceleracion. El cambio fundamental 
es el de transformar el concepto de fuerza de atraccion al de deformacion (o mejor curvatura) del 
espacio-tiempo, y asf una particula no experimenta una aceleracion por una fuerza sino una trayectoria 
en geodesicas en un espacio-tiempo curvo. El ente que realiza esta curvatura en el espacio-tiempo es 
cualquier forma de materia-energia, de modo que siguiendo [23]: La curvatura le dice a la materia como 
moverse, pero la materia le dice al espacio como curvarse. 

2.2. Postulados y Formulacion Matematica de la Relatividad General 

Daremos entonces un repaso de los elementos basicos de la RG siguiendo [23]- [30]. Si bien Einstein 
logro encontrar una teoria geometrica de la gravitacion, su construccion se basa en postulados, que sin 
embargo fueron desarrollados buscando que las ecuaciones dinamicas den los Hmites correctos (gravedad 
Newtoniana y Relatividad Especial). El primer postulado se basa en la descripcion del espacio-tiempo: 

Postulado 1: El espacio-tiempo esta descrito por un par {A4, g) siendo A4 una variedad ^-dimensional 
y g una metrica Lorentziana sohre A4. 

La curvatura en la variedad esta descrita por el tensor de curvatura de Riemann R, que en componentes 
^La equivalencia entra masa inercial y masa gravitacional de un cuerpo. 
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Capi'tulo 2. Elementos de Relatividad General 



puede escribirse como 

= ^T^'sp ~ 95^"l3 + '^yrr'^SlS ~ ^^^r^^, (2.1) 

con las conexiones (o simbolos de Christoffel) 

= ^5"" [d^gap + dpg^^ - d^g^^] , (2.2) 

en donde gafi son las componentes del tensor metrico. Una restriccion fundamental que se impone sobre 
la variedad M en RG es que sea libre de torsion, lo cual se basa en el siguiente teorema [31] 

Teorema 2.2.1 (Teorema Fundamental de la Geometria Riemanniana) . En una variedad Riemanniana 
{A4,g), existe una unica conexion simetrica F^^ que es compatible con la metrica g. Esta conexion es 
llamada la conexion de Levi-Civita. 

La condicion de que la variedad sea libre de torsion se satisface exigiendo que el tensor de torsion de 
Cartan S^^ sea nulo 

^p-y = r^7] = 0' (2-3) 

y la relacion entre la conexion F^^ y la metrica g se obtiene imponiendo que la derivada covariante de 
la metrica sea nula 

V75a/3 = 0, (2.4) 
con la definicion usual para la derivada covariante 

El postulado de conservacion de energia juega un papel fundamental en RG. 

Postulado 2: Conservacion local de la energia: Existe un tensor simetrica T^p = Tap{'ip) = Tj^^ 
que es funcion de los campos de materia ^ y sus derivadas tal que : 

I. Tq,^ = sobre U C Ai si y solo si ipi = para todo i sobre U. 
II. V^T"^ = 0. 

Finalmente, se encuentra el postulado de las ecuaciones de campo de Einstein, las cuales nos van a dar 
las relaciones entre la geometria (curvatura) con los campos de materia: 

Postulado 3: La metrica sobre la variedad espacio-tiempo (A^,g) esta determinada por las ecuaciones 
de campo de Einstein 

(2.6) 

siendo Rai3 el tensor de Ricci (Rai3 = Rarns)' ^ escalar de curvatura (R = g'^^Rap), T^p el tensor 
de energia-momentum, y n = SvrG, G la constante de gravitacion universal y usamos unidades de c= 1. 



per rpOll, 



(2.5) 



Raf5 — 7,R9a(5 — K-Ta 



2.2. Postulados y Formulacion Matematica de la Relatividad General 
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Si definimos el tensor de Einstein Gap como 

Gap = Rap — ^Rga/^, (2-7) 

entonces se tiene una restriccion geometrica 

V^G"^ = 0. (2.8) 
ecuacion conocida como Identidad de Bianchi. 

El movimiento de una particula es descrito por su trayectoria en el espacio-tiempo, x"(A), donde A 
es un parametro. Una particula libre, i.e., una particula sobre la que no se ejerce fuerza alguna (otra 
diferente de la gravedad), satisface la ecuacion de la geodesica, la cual se escribe 

y^V^y^ = 0, (2.9) 

donde = es el vector tangente a la trayectoria. Esta ecuacion puede escribirse en la forma 
conocida 

La expresion aplica para los siguientes casos: 

■ Particulas con masa, en cuyo caso usualmente se toma como parametro A el llamado tiempo 
propio tal que el correspondiente vector tangente esta normalizado: ga^y^V^ = — 1- 

■ Particulas sin masa, en particular el foton, en cuyo caso el vector tangente, usualmente denotado 
por /c" es nulo, i.e., ga/sk'^k^ = 0. 

Las ecuaciones de campo de Einstein nos determinan entonces la metrica dado un contenido de materia- 
energfa, y las ecuaciones del potencial gravitacional Newtoniano (la ecuacion de Poisson) se encuentra 
como un Ifmite de tales ecuaciones. La relatividad especial es entonces un caso particular de la relatividad 
general, para la cual M es una variedad conformalmente plana y g = t/ con rj el tensor de Minkowski. 



2.2.1. Principio Variacional en Relatividad General 

Un problema que enfrento Einstein despues de escribir su Teoria de la Relatividad General fue el de 
encontrar sus ecuaciones de campo a partir de un lagrangiano y una principio variacional 6S = con S 
expresando la accion total. Einstein y Hilbert encontraron, de manera independiente, la accion asociada 
al campo gravitacional, la cual es denominada como la accion de Einstein-Hilbert. En terminos esta 
Seh, el termino de frontera de Gibbons- York-Hawking Sgyh (el cual fue introducido para relajar las 
condiciones de frontera como veremos mas adelante) [32], [33], y la accion asociada con todos los campos 
de materia 5a/, la accion total se puede escribir como [25]: 

S = :^{Seh + Sgyh) + Sm, (2.11) 
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donde 



Seh = 
Sgyh = 2 



(fx 



-gR, 



d^ye 



dV 



(2.12) 



(2.13) 



aqui V es un hipervolumen en A4, dV su frontera, h el determinante de la metrica inducida, K es la 
traza de la curvatura extrmseca sobre la frontera dV, y e es igual a +1 si dV is como de tiempo y — 1 if 
dV es como es espacio (se asume que dV no es nulo en alguna parte). Las coordenadas son usadas 
para la region finita V y para la frontera dV. Vamos a obtener ahora las ecuaciones de campo de 
Einstein variando la accion con respecto a g"'^. Fijamos tal variacion a la condicion [24], [25] 



0, 



(2.14) 



dV 



es decir, la variacion de la metrica se anula en la frontera dV. Usamos ahora los siguientes resultados 
[25], [28] 

Sgap = -ga^gpJg"", 5^"^ = -g'^'g^'Sg^u, 
= -^V^gafiSg"'^, 
v^(5r^^)-V5(5r-^), 



V/3(5n 



(2.15) 

(2.16) 
(2.17) 
(2.18) 



(2.19) 



(2.20) 



/3a' " P\"'' -yaJ- 

Daremos a continuacion un detallado repaso de los principios variacionales en RG siguiendo para esto 
[24], [25]. La variacion del termino de Einstein-Hilbert es 

SSeh= I d'^x{R6^/^ + ./^6R). 
Jv 

Ahora con R = g'^^Ra/B, tendremos que la variacion del escalar de Ricci es 

6R = 6g''^Rap + g''^6Rap. 

usando la identidad de Palatini (2.18) podemos escribir [28] 

6R = <55"^i?a/3 + 5"^(V^(<5r^J - Vp{5n^)), 

= 6g"^Ra(s + V45"^(5r^J - g'^^'iSTl^)), (2.21) 

donde hemos usado la compatibilidad metrica V^ga/s = y renombramos algunos indices mudos. 
Reemplazando estos resultados para las variaciones en la expresion (2.19) tenemos 

SSeh= [ d^x{R5^/^ + ^/^5R), 
Jv 

= ^ d'^x (^-^Rg^^^Sg-^ + Ra^V^6g"^ + ^g^^ {SV^^^) - ^"'^(^r^^))) , 

= j^d^x^gl^Rap - \Rga^5g"P + j^d^x^gV,{g^P{6Tl^)-g'''^{5Vl^)). (2.22) 
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Denotando el termino de divergencia por 6Sb, 

6Sb= I d''x^V^{g^^{5T'^p,)-g"^{5Tl^)), 
Jv 



definimos 



de modo que el termino de frontera se puede escribir como 

SSb= [ (fx^gV„V\ 
Jv 

Usando el teorema de Gauss-Stokes [25], [28] 

Jv JdV 



(2.23) 
(2.24) 

(2.25) 
(2.26) 



donde es un vector normal unitario a dV. Usando esto podemos escribir (2.25) como el siguiente 
termino de frontera 

5Sb= i d^ye^/WWaV, (2.27) 

JdV 

con V'^ dado en (2.24). La variacion (^F^^ se obtiene usando que F^^ es el simbolo de Christoffel {^q,} 



(2.28) 



obteniendo 



= \^9'''^ [dp9ia + dag^p - d^gpa] + ^g"'^ [d^{6g^a) + da{5g^p) - a^(%a)] . 
De la condiciones de frontera Sgap = 5g°^^ = la variacion (2.29) es 



_ 1 

dv~ 2 



g"'^ [di3{5g^a) + da{6g^(s) - d-y{6gi3a)] ■ 



podemos escribir 



dV 



'-^g^"'[dp{5g^a) + da{5g^p) - 5^(55/3a)] 



1^7 ) 



Vn 



av 



dV 



^g^'^ [dp{5g-ia) + da{5g-fp) - d^{5gpa)] 
g'^'^daiSg^y 



(2.29) 

(2.30) 
(2.31) 

(2.32) 



^djg^^ [dfsidg^a) + da{6gy^) - d,{Sgpa)] - Is^g"^ [daiSg 
g''^[dp{dg,a)-d^{6gp^)]. 



(2.33) 
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Ahora evaluamos el termino n'^Va]^^ usando para esto que 

= /i"^ + en"n^, (2.34) 

entonces 

n-y, = n^h''^ + en''n'')[dfs{5g^a) - 5.(%„)], 
aV 

= rfh''^[dfi{5g^^) - d^{5gpo)], (2.35) 

donde usamos la parte antisimetrica de en°'n^ con s = n'^n^ = ±1. Ahora, sabiendo que Sga/s = en 
la frontera tendremos que /i"'^5^((5g(Ta) = [25]. Obtenemos finalmente 



" av 



-n''h^^d,{6gpa)- (2.36) 
De modo que la variacion del termino de Einstein-Hilbert es 

5Seh = 2 j^<fx^g{Ro^p - ^i^ffa/?)*^^"^ - £^d^y ey^\h''''d^{Sg^aK . (2.37) 

Como podemos observar, si fijamos Sgap = existe un termino de frontera adicional. En la mayoria de 
textos tal termino se anula argumentando flujos nulos en el infinito; otro argumento es el de fijar tanto 
la variacion de la metrica como su primer derivada a ser nulas sobre la frontera, es decir, 6gai3 = y 
d-ySgai3 = 0. Aunque este ultimo argumento nos conduce directamente a las ecuaciones de campo de 
Einstein (pues la contribucion en la frontera es automaticamente cero), implica fijar dos condiciones 
en la variacion. Para evitar esto Hawking, York y Gibbons introdujeron un termino de frontera que 
permite tener un problema variacional bien definido con solo 5gai3 = 0. Es importante recalcar que 
aunque la accion total se vea modificada por terminos adicionales sobre la frontera, las ecuaciones de 
campo obtenidas son las mismas. 

Consideramos entonces la variacion del termino de frontera de Gibbons- York-Hawking 

6Sgyh= i d^ye^M^K. (2.38) 
Jav 

Usando la definicion de la traza de la curvatura extrfnseca [25] 

K = Van", 

= (/i°^ + en''n^)Vf3na, 
= h'^Npna, 

= h"^{dpna-T}^n^), (2.39) 
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la variacion es 

= \h^^d^{5gp^)n\ (2.40) 

Esto viene de la variacion (5r^^ evaluada en al frontera, y del hecho que h°^^ dp{5gtjoi) = 0, h'^^da{Sgai3) = 
0. Entonces tenemos para la variacion del termino de frontera 



5SGYH=i d'yey/\h\h''^d„{59Pc.)rf. (2.41) 
JdV 



Vemos que este termino cancela exactamente la contribucion en la frontera proveniente de la accion de 
Einstein- Hilbert. Ahora, si tenemos una accion de materia definida por 



Sm = d x^/^CM[ga|3,'4^], (2.42) 
Jv 

donde '0 denota todos los campos de materia. La variacion de esta accion toma la forma 

6Sm= / d^x6{^/^CM), 
Jv 



d'xl^6g'^f'V^ + CM6V^ 



d^^i-Q^- 2^M9a^ I'Js"^ (2.43) 



como es usual, definimos el tensor de energia-momentum por 



_ ^dCM ^ 2 6Sm ... 

entonces: 

^SM = -\j^ d^x ^gT^pSg''^, (2.45) 

imponiendo que las variaciones totales permanezcan invariantes con respecto a 5g°^^ , podemos escribir 
finalmente 

\ S S 1 

0, =^ Ral3 - -RgajS = I^Ta/B, (2.46) 



que corresponden a las ecuaciones de campo de Einstein (2.6). 
2.2.2. Relatividad General a la Palatini 

Un metodo alternativo para encontrar las ecuaciones de campo fue introducido por Palatini en 1919. 
Su formulacion consiste en tratar la metrica y la conexion T'^^ como dos campos independientes. 
Para esto definimos entonces un tensor de Ricci TZa0 definido por 

= Kr,^ = d.rlp - dpT^va + - Kfs^va^ (2-47) 
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y el escalar de Ricci como TZ = g"TZai3- La accion de Palatini puede escribirse como 



Sp 



1 

2k 



(2.48) 



con Sm siendo nuevamente la accion asociada a los campos de materia (2.42). Usando los resultados 
(2.16) y (2.21) asumiendo que la conexion es independiente podemos escribir para la variacion con 
respecto a gaS 



a/3 



(2.49) 



e imponiendo que la accion sea invariante con respecto a 5gap obtenemos que 

1 55p ^ _ 1 



(2.50) 



que corresponden a las ecuaciones de campo de Einstein. Consideramos ahora la variacion con respecto 
a usando para esto la identidad de Palatini (2.18) 



5Sp = ^ f d'x^gg^P{V,{5TlJ - V p{5Tl^^)), 

^K Jy 



(2.51) 



integrando por partes, y tomando en cuenta que <^r^^ = en al frontera, de modo que terminos lineales 
en (^r^^ son cero, podemos escribir 



^^P = -h I d'x{V,{^gg'^P)5Tlp-Vp{^gg^f')5^:), 
55p = -^ / (V7(^/=59"^)-V,(V^5"'^)5^)5^^^, 
imponiendo de nuevo que esta variacion sea invariante con respecto a tendremos 



(2.52) 



tomando la traza de esta expresion tenemos 



-gg 



0, 



(2.53) 



(2.54) 



que es la condicion de que la conexion sea de Levi-Civita. Con el metodo de Palatini se recuperan 
entonces las ecuaciones de campo y ademas se obtiene la relacion entre la metrica y la conexion. 



Las ecuaciones de campo de Einstein (2.6) son entonces el punto de partida para la evolucion dinamica 
de las cantidades fisicas en un espacio-tiempo. 



2.3. Efectos de Curvatura = Ecuacion de Desvfo Geodesico 
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Figura 2.1: Desvi'o Geodesico. 



2.3. Efectos de Curvatura = Ecuacion de Desvfo Geodesico 

Una forma elegante de observar los efectos de la curvatura en una variedad es a traves de la denominada 
Ecuacion de Desvw Geodesico [34]. Vamos a continuacion a hacer una descripcion de tal ecuacion 
siguiendo [24]-[26]. Sean 70 y 71 dos geodesicas vecinas con parametro afm i^. Introducimos entre las 
geodesicas una familia de geodesicas que se interpolan con parametro afm s, Figura 2.1. El campo 
vectorial V" = es tangente a las geodesicas con parametro affn z^, y la familia s tiene a ry" = 
como su campo vectorial tangente. Vamos a derivar la expresion para la aceleracion del vector desviacion 
7]°'. Para esto consideremos las coordenadas Xq para la geodesica 70 y para 71. La evolucion de la 
geodesica 70 base viene descrita por la ecuacion 

y scan x" = Xq + 77° las coordenadas de la geodesica vecina, la cual satisface 

d^xf dx[dxl_ 

du^ dv du ^^-^^^ 

La conexion afm a primer orden se puede escribir 

T%{xt) = T%{x^ + r?") - r^,(x?) + d^T'l^ix^W, (2.57) 

por otro lado 

dx? dx?. dri" 
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de modo, la ecuacion resultante para la ecuacion (2.56) es 

Reescribiendo Xq = x", ^'^^{xq) = FJ^^, tenemos 

d^Tj- d^x- (dx^ dr^P\ (dx-^ drj-y 



(2.60) 



d^ry" / d'^x" dx^ dx'^ \ / dx^ dr]'^ dx"' dr]^ 

~d^^\~d^^ + V ^^^d^^ ^^^^ 



El termino del primer parentesis se anula debido a la ecuacion geodesica (2.55). Los dos terminos del 
segundo parentesis se pueden escribir en un solo termino, debido a la simetria de la conexion, de modo 
que 

^ + a r" 77'^— — + 2r- -0 f2 62) 

Esta ecuacion es posible escribirla de manera mas compacta acudiendo al concepto de derivada cova- 
riante a lo largo de una curva. Sea un campo vectorial definido sobre una curva cuyo parametro 
afin es v. Al igual que como se define la derivada covariante sobre todo el espacio, la operacion definida 
como 

DA' ''^%r?,^^., (2.«3) 



Du du du 

du 



define la derivada covariante a lo largo de la curva, con Vl^ = el vector tangente a la curva. Escrito 
de forma mas compacta 

= V^A'^V^. (2.64) 
Como ^^jj^ es otro campo vectorial, podemos tomar su derivada covariante a lo largo de esta curva 

utilizando la definicion de derivada covariante tenemos 



entonces 

[-D^) = ^ KsV-A') + r,",^^^- (2.67) 

Por otro lado tenemos que 

^ /dA"\ d fdA"\ dx"/ d fdA''\ d^A" 
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de modo que 



7/3 



dAl^ 
du 



VlpV'AP^ V^. (2.69) 



Debido a que V"' es el vector tangente a la curva geodesica, esto es, la derivada covariante de a lo 
largo de la misma curva del vector tangente es nula, = 0, podemos utilizar 

d^VV^ = -Tl^VV^, (2.70) 



de esta forma 



■ aS u-y 



Ahora, volviendo a la ecuacion (2.62), ecuacion para la separacion de dos geodesicas proximas, despe- 



jando ^ 



du^ 



du^ ' 



(2.72) 



De esta manera, si consideramos a r?" como el campo vectorial A" , reemplazando el valor de la derivada 
parcial segunda de la ecuacion (2.71) y reescribiendo algunos indices mudos, tenemos 



n2 a ^ p 



dr]^ 



du 



(37 •'•''/ fj^'- up 

_i_ \n yS _|_ -pa 



du 



+ rf^yV)^^ (2.73) 



agrupando terminos semej antes 



Du^ 



+ 



dr]^ 



,dri^ 



djf 



T^VlpT^^VVP + r°«rf,r?^l/'^y^ (2.74) 



dado que F^^ = V^g (torsion nula), los terminos del segundo parentesis se anulan y 



-7/3 



Du^ 



(2.75) 



2^a 



-dsT^r^^VV^ + dsT^^r^^'vW^ - C^r^^T]^^^^ + T^^r^^r/^y (2.76) 



donde de nuevo renombramos algunos indices mudos. Cambiando en el primer termino despues de la 
igualdad 5 < — > 7, en el tercer termino 5 < — > 7, y en el cuarto termino 5 < — > 7 y /3 < — > 7, tendremos 



(2.77) 
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El termino entre parentesis es el tensor de curvatura de Riemann (2.1), por lo que la ecuacion de desvio 
geodesico queda escrita como 

(2.78) 

Esta importante ecuacion nos muestra que la aceleracion del vector desviacion depende del tensor de 
Riemann, y por lo tanto da una directa prueba de la curvatura en una variedad. Como es de esperarse 
para un espacio-tiempo piano R'fj^s = 0, y por lo tanto el vector desviacion no sufre aceleracion. Esta 
importante relacion nos permitira describir todas las propiedades geometricas de un espacio-tiempo, y 
nos dara de una forma elegante las cantidades observacionales (en particular las distancias) [35]. 



2„a 
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Capitulo 



3 



El Modelo Estandar de la Cosmologia 



a Relatividad General (RG) es la base para la construccion de modelos que describan la evolucion 



del universe. Despues de ser verificada experimentalmente, se empezaron a construir modelos a 
partir de las ecuaciones de campo de Einstein. Una vez dado el fracaso del modelo estatico de Einstein 
motivado por observaciones en el corrimiento al rojo de galaxias, surgieron nuevos modelos basados 
en las propiedades estadisticas del uni verso a grandes escalas. El acunado actualmente como "Modelo 
Estandar de la Cosmologia" reposa sobre dos principios fundamentales: la homogeneidad e isotropia. 

Estc universe liomogeneo e isotropico descrito por la metrica dc Robertson- Walker y denominado co- 
mo cl univcrso dc Fricdamnn-Lcmaitrc-Robcrtson- Walker (FLRW) cs la base actual dc la cosmologia 
modcrna. Aunquc cs accptado y varias prucbas cxpcrimcntalcs lo soportan, cxistcn varios modelos que 
guardan alguna de las propiedades del denominado principio cosmologico, entre ellos los universes de 
Bianchi (homogeneos pero no isotropicos) y los de Lemaitre-Tolman-Bondi (isotropicos pero no homo- 
genos). 

En lo que sigue daremos los elementos basicos del modelo estandar de la cosmologia centrandonos 
en sus problemas actuales y en el papel fundamental de la medicion de distancias. 

3.1. La metrica de Robertson-Walker y la expansion del universo 

Como mencionamos anteriormente, el Modelo Estandar dc la Cosmologia esta basado en dos principios: 
La homogeneidad e isotropia del universo. Estos postulados permiten escribir el elemento de li'nea en 



"An observer situated in a nebula and moving with the nebula 
will observe the same properties of the universe as any other 
similarly situated observer at any time. " 

Sir Hermann Bondi, 1948. 
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un universo maximalmente simetrico por medio de la metrica de Robertson- Walker [2] , [36] , la cual esta 
dada por 



dr^ 



1 — kr"^ 



+ r'^de^ + sin^ edip"^ 



(3.1) 



siendo a{t) el factor de escala y A; un parametro que nos dice si la geometria del universo es plana 
{k = 0), esferica (A; = 1) o hiperbolica (k = —1), Figura 3.1. La variedad ahora esta representada 
por una foliacion = S x R en donde S representa una tri-variedad maximalmente simetrica y M la 
coordenada temporal del espacio-tiempo. El principio cosmologico permite ademas escribir el tensor de 
energia-momentum en la forma de un fluido perfecto 



Tap = (P + P)UaUf3 + pg^/^, 



(3.2) 



siendo p la presion del fluido, p la densidad de energia y la cuadrivelocidad de los observadores 
fundamentales. Las ecuaciones de campo de Einstein (sin constante cosmologica) para la metrica (3.1) 
con el tensor de energia-momentum (3.2) son 



K k 
3 



Componente temporal 



(3.3) 



a 
a 



(p + 3p), 



Componentes espaciales 



(3.4) 
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en donde definimos el parametro de Hubble como H = ^. ha conservacion de energia, expresada por la 
ecuacion V^r°^ = se escribe para esta metrica como 

Para resolver estas ecuaciones es necesario conocer una relacion entre la presion p y la densidad de 
energia p. La forma mas usual es considerar una ecuacion de estado de la forma [36] 

p = wp, (3-6) 

con w un parametro, que a lo mas puede depender del tiempo w = 'w{t). Al conjunto de ecuaciones 
(3.3)-(3.5) se les conoce como ecuaciones de Friedmann-Lemaitre, y junto con la ecuacion de estado 
(3.6) determinan el modelo cosmologico. Los modelos mas trabajados son los de polvo (ti; = 0) y 
radiacion {w = 1/3), los cuales dan cuenta de la evolucion del universo. Sin embargo, hoy en di'a 
sabemos que nuestro universo se expande de forma acelerada. De esta manera surgen los modelos de 
constante cosmologica, que tratan de explicar esa expansion introduciendo formas exoticas de energia. 
Para que tengamos un universo en expansion se necesita que a > 0, y a partir de la ecuacion (3.4), 
junto con la ecuacion de estado (3.6), se llega a la condicion w < —1/3. Los modelos con ecuaciones de 
estado de esta forma se denominan de Quintaesencia [4]. El caso con w = —1 se obtiene suponiendo en 
la ecuacion (3.5) que la densidad es constante, lo que nos lleva a una ecuacion de estado de la forma 

P = -P, (3.7) 

es decir, una ecuacion de estado con presion negativa. A esta forma de energia que reproduce una presion 
negativa se le conoce como energia oscura. Podemos entonces asociar una densidad, que llamaremos 
densidad de energia del vacfo relacionada con la constante cosmologica por la relacion [9] 

PA = -, (3-8) 

K 

que se puede obtener si introducimos el termino asociado a la constante cosmologica en la ecuacion 
(2.6) ^ 

Rap - 2^90/3 = KTai3 + Ag^p- (3.9) 
Estas ecuaciones de campo nos conducen a la siguiente ecuacion de Friedmann 

H^ = -p-^ + -, 3.10 
3 3 

definimos los parametros cosmologicos por las siguientes relaciones 

^m^^, ^^A=3^, ^'^-Jh^^ (^-^^^ 

con lo cual la ecuacion primera ecuacion de Friedmann se escribe como 

f)^ + J7A + J]fc = 1, (3.12) 



3.1. La metrica de Robertson-Walker y la expansion del universo 
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No Big Bang ^ ..,99% 




12 3 

Figura 3.2: Regiones de coincidencia para los parametros JIa y H™, (Perlmutter et al, (1999)) [3]. 



llamada como la regla de suma cosmica. Observaciones con supernovas de tipo la (Sn la) muestran que 
los valores actuales de estos parametros son [3], [37], [38] 

K. 0.28, VLk ^ 0.72, nk ^ 0. (3.13) 

La figura 3.2 muestra las regiones de coincidencia en un piano de y a partir del cual se pueden 
estimar las cotas para los anteriores parametros. La contribucion al parametro de curvatura 0,^ es nulo 
y + = 1- Por otro lado, gracias al estudio de la Radiacion Cosmica de Fondo (CMB) por medio 
del proyecto Wilkinson Microwave Anisotropy Probe (WMAP), se demuestra que el contenido actual 
de materia en el universo esta compuesto en un 72.1 % de energfa oscura, un 23.3% de materia oscura 
y un restante 4.6% de materia ordinaria (o mejor llamada barionica)^ [37]. 

Asi pues nos encontramos actualmente en un universo que se expande aceleradamente y con abun- 
dancia de un contenido de materia y energfa que aun desconocemos. Aunque hoy en dfa se maneja el 
modelo ACDM (constante cosmologica + materia oscura frfa), se tienen sin embargo indicios que nos 

^Obtenidos tambien de http://map.gsfc.nasa.gov/. 
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llevan a reconsiderar la introduccion de la constante cosmologica. Dados los valores para los parametros 
cosmologicos se puede dar una aproximacion, al menos en ordenes de magnitud, para la densidad de 
energia relacionada del orden de 2 x 10"^'' erg/cm^. Sin embargo a nivel de teorias cuanticas, la energia 
del punto cero de los campos, la cual se puede asociar a la constante cosmologica, es del orden de 
2 X 10^^*^ erg/cm^. Esta diferencia de 120 ordenes de magnitud para los valores experimentales y teori- 
cos para la constante cosmologica es el conocido problema de la constante cosmologica. Otro problema 
fundamental es el denominado problema de coincidencia, y es el por que actualmente la densidad de 
energia de la constante cosmologica es comparable a la de la materia [9] . 

Es por estos problemas y la imposibilidad de dar una explicacion razonable a la energia oscura que 
se ban propuesto teorias alternativas para explicar la actual expansion acelerada del universe. Una 
de ellas es considerar teorias de gravedad de alto orden, en las cuales se pueda ver tal expansion co- 
mo resultado de terminos adicionales en las ecuaciones de campo de Einstein, debidos solamente a 
contribuciones geometricas. 



3.2. Distancias en el Modelo Estandar 



Un aspecto fundamental en el modelo estandar de la cosmologfa se refiere a la medicion de distancias. 
Como veremos tal medicion depende de la metrica usada y tambien del contenido de materia en el uni- 
verso. Como vimos anteriormente, las ecuaciones de Friedmann dependen principalmente de la forma 
de las ecuaciones de campo, de modo que las expresiones para las distancias dependeran finalmente de 
la forma de la metrica y las ecuaciones de campo usadas. 

Las observaciones astronomicas estan basadas en la radiacion que viaja a traves del espacio-tiempo 
y que llegan a nosotros a traves de geodesicas nulas. En el caso del universo de FLRW es posible 
escoger geodesicas radiales, por lo que para estas geodesicas ds^ = d9 = dcj) = 0, y por lo tanto 
= —dt'^ + a'^{t)dr'^. De aqui se sigue que la radiacion emitida en un punto E y recibida en otro punto 
O obedece las siguientes relaciones 

O f^o dt ro da 

con a{t) = 

3.2.1. Redshift 

Vamos a considerar un pulso de luz que se emite en n punto E y llega a un punto O. Definimos 1 + Zc, 
siendo Zc el redshift cosmologico, como la razon entre la longitud de onda medida por el observador O, 
\o y la longitud de onda emitida en A^;. Dado que la longitud de onda es inversamente proporcional 
al pen'odo entonces 

(' + ") = A^ = AT^- 
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Podemos interpretar el periodo AT como el tiempo entre cresta y cresta de la onda de luz. Ahora los 
periodos ATq y ATe se relacionan con el factor de escala por [39], [40] 

ATo ATe 



a{to) a{tE) ' 

de esta manera 



(3.16) 



(1 + ^) = ^ = ^, (3.17) 
a{tE) 

en donde a(to) es el tamano del universo en el tiempo en el que la luz del objeto es observada, y a{te) 
es el tamano en el tiempo en que la luz fue emitida. 

Si consideramos que a (to) corresponde al factor de escala evaluado actualmente de modo que pode- 
mos fijar a(to) = 1 y con a{t) el redshift evaluado en algiin tiempo cosmico correspondiente a un 
redshift z, la expresion anterior se reduce a la relacion conocida 



a{t) = (3.18) 



3.2.2. Distancia Comovil x{z) 

Para escribir las distancias de forma general vamos primero a escribir la ecuacion de Friedmann (3.10) 
considerando la densidad p compuesta por una contribucion de materia pm y otra de radiacion p^ 



H' = -iPm + Pr)-^ + -, (3.19) 

3 3 

estas contribuciones satisfacen las ecuaciones de estado pm = y pr = gPr, de modo que las densidades 
de materia y radiacion satisfacen 

Pm oc a~^, pr oc a~^, (3.20) 

relaciones que se deducen a partir de la ecuacion de conservacion (3.5). Podemos escribir entonces el 
parametro de Hubble de la siguiente forma 

H\a) = Hl[Vlmoa-^ + ^roa'^ + S^fcoa"^ + ^^a), (3.21) 

siendo Hq el valor actual para el parametro de Hubble y definimos ademas las siguientes cantidades 
(ver ecuacion (3.11)) 

Podemos escribir el parametro de Hubble en funcion del redshift usando (3.18) 

H^{z) = Hi (J7„o(l + zf + VLro{l + zf + 17fco(l + z)a^ + J^a) , (3.23) 
con estas cantidades definimos la distancia comovil xiz) como 

(3.24, 

Esta distancia es la mas fundamental en cosmologia, pues como veremos las otras relaciones estan 
derivadas en funcion de esta. 
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Figura 3.3: Distancia diametral angular. 



3.2.3. Distancia Comovil (transversal) dM{z) 



Se define la distancia comovil transversa dM por las relaciones 

1 



HQ\/\yLk 



smh[^/^^Hoxiz)], ilk > 0, 



S^fc = 0, 



siny-^kHoxiz)], < 0, 



(3.25) 



esta distancia tambien es conocida como distancia de movimiento propio definida como la razon entre 
la velocidad transversal y el movimiento propio [1]. 

3.2.4. Distancia Diametral Angular d^iz) 

La distancia diametral angular dA es definida como la razon del diametro propio de un objeto D con 
su diametro angular aparente 6 (en radianes), Figura 3.3. 

D 



dA 



(3.26) 



de la metrica de Robertson- Walker (r = cte y (p = cte), el diametro propio de la fuente con A6 = 5 es 

D = a{ti)ri5, (3.27) 

de modo que la distancia diametral angular queda 

dA = a{ti)ri. (3.28) 
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Esta distancia esta relacionada con la distancia comovil transversa por 



dA 



(3.29) 



de modo que usando la expresion (3.25) tendremos 

r 1 



dAiz) 



1 + z 



Ho\/Wk 



smh[^/^^Hox{z)], > 0, 



: siny-nkHoxiz)], < 0, 



(3.30) 



3.2.5. Distancia de Luminosidad dL{z) 

Finalmente la distancia de luminosidad esta definida por la relacion entre el flujo de una fuente 
y su luminosidad L 

d, = yx. (3.31) 

La relacion de la distancia de luminosidad con la distancia comovil transversal y la distancia diametral 
angular por 

dL = {l + z)dM = (1 + z)^dA, (3.32) 

de modo que 

r 1 

sinh[vn^i^ox(^)], ^k > 0, 



dLiz) = (1 + z) < 



Xiz), 



Qfe = 0, 

siny-QkHoxiz)], 0.k < 0, 



(3.33) 



Las relaciones anteriores para las distancias nos permiten comparar entonces las observaciones con los 
modelos teoricos, y en especial determinar los valores de los parametros cosmologicos. 



3.3. Ecuacion de Desvio Geodesico en FLRW 



Una forma muy elegante de tratar el problema de distancias en cosmologia (y en general todas las 
propiedades del espacio-tiempo), es a traves de la EDG. Daremos entonces a continuacion un repaso de 
la EDG en el universo FLRW siguiendo [35], [39]. 
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El punto fundamental es relacionar las cantidades geometricas de las ecuaciones de campo (el ten- 
sor de Rice Rap y el escalar de curvatura R) con el tensor de energia- momentum Tq,^. Para esto vamos 
a considerar la siguiente descomposicion del tensor de Riemann en sus partes libres de traza y con traza 
[24], [41] 

1 R 

Ra/B'jS = Ca(3y5 + ^ {da-yRs/S — daSRyfi + dpsR-ya — dp-yRsa) — {gaydSp — 9a&9yl3) , (3.34) 

en donde Cap-yS es el denominado tensor de Weyl. A partir del elemento de Imea (3.1) se muestra que 
el tensor de Weyl es nulo para universos homogeneos e isotropicos CapyS = ^. A partir del tensor de 
energia-momentum (3.2) podemos obtener su traza 

T = 3p- p. (3.35) 

y de las ecuaciones de campo de Einstein (con constante cosmologica) 

Ral3 — -^Rgap + ^9ap = I^Tafj, (3.36) 

podemos escribir el escalar de curvatura -R y el tensor de Ricci R^p como 

R = k{p- 3p) + 4A, (3.37) 

Rap = K{p + p)UaUp + ^[K{p-p) + 2A]gap, (3.38) 

de modo que reemplazando estos resultados en (3.34) obtenemos 

RaPyS = ^k{p + p)[gayUpUS - QaSUpU^ + QpsUaU^ - Qp^UaUs] 

+ ^(/ip + h)[ga^gp8 - gasgpy]- (3.39) 

Si el campo vectorial V" esta normalizado, VaV" = e, tenemos 

Rap^sV^V^ = ^K{p + p)[gayUpusV^V^ - gasV^upu^V^] 

+ ^K{p + p)[gpsV^V^UaU^ - gpyV^V^UaUs] + \{^p + K)[gaygp5V^V^ - gaSgpyV^V^l 
de donde obtenemos finalmente que 

RaPysVf^V^ = \{Kp + K)[egay-VaV^] + ]^k{p + p)[{upV^fga^ - 2{Vpuf^)u^aVy) + eu^u^l (3.40) 
subiendo el primer mdice del tensor de Riemann 
RlysV^V = hnp + K){e5^.^ - V^V,) + \n{p + p)[{Vpu^)H"^ - {Vpu^)u^V, - {Vpu^Wu,], (3.41) 



^En universos con anisotropfas, como por ejemplo los de Bianchi, el tensor de Weyl es en general no nulo. 
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multiplicando por r]'^ 



- {Vpu^)u"{v^l^'') - {v^u^)v-{u^l^'')] (3.42) 

usando ahora que "qau" = y r]aV'^ = obtenemos finalmente [35] 



7] 



(3.43) 



donde E = — V^u". Esta es conocida en la literatura como la ecuacion de Pirani [42]. Notemos que este 
termino fuerza es proporcional al mismo r/", i.e., de acuerdo con la EDG solo la magnitud rj cambia 
a lo largo de la geodesica, mientras que su orientacion espacial permanece fija. Consecuentemente la 
EDG (2.78) se reducira simplemente a una relacion diferencial para la cantidad escalar rj. Esto refleja la 
isotropia espacial del tensor de curvatura de Riemann alrededor de cada punto en la presente situacion. 

Consideraremos a continuacion dos casos particulares de la GDE que nos permitiran determinar tambien 
las distancias cosmologicas en los universos de FLRW. 



3.3.1. EDG para Observadores Fundamentales 

En este caso, el campo vectorial geodesico es el campo vectorial de la cuadrivelocidad del fluido m°. 
El parametro affn coincide con el tiempo propio de observadores fundamentales, i.e., ly = t. Como estas 
son geodesicas temporales, e = — 1 y dado que los campos vectoriales estan normalizados, £^ = 1 de la 
ecuacion (3.43) se tiene 



^(p + 3p) + ^A 



Si el vector desviacion es r]a = iea, la isotropia implica que 



Dt 



0, 



tenemos ademas que eQ,e° = 1, ecU^ 



0. Con base en esto, segiin (3.44) y la EDG 



1^ , 1 . 

-{p + 3p)--A 



Dado que la derivada covariante cumple la regla de Leibnitz 



Dt]"' _ D 
1m ~ Dt'' 

D'^ri" 



D£ ( De"" 

Dt V Dt 



D£ 
Dt' 



Dt V Dt 
Dh 



D' 



Dt'^ 



DiDe"^ 



Dt'^ Dt"^ 



(3.44) 



(3.45) 



(3.46) 



(3.47) 



(3.48) 
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de modo que 



^(p + 3p)-^A 



(3.49) 



Esta es la ecuacion de Raychaudhuri [35]. Escogiendo entonces £ = a, siendo a el factor de escala de 
modo que la ecuacion de Raychaudhuri se reduce a 



a K , s A 
a D 3 



que corresponde a la ecuacion de Friedmann (3.10). Multiplicamos por aa 

dd H — k{p + 3p)ad Aad = 0, 

6 3 

y usando la ecuacion de conservacion de la densidad de materia (3.5) podemos escribir ^ 



-3-k{p + p), 
a 



de esta manera, al considerar la derivada de Kpa^ 

d^npa?') 



dt 



a np + 2anpd, 



-3-K,{p + p) I +2aKpd, 



reorganizando terminos 
Reemplazando (3.54) en (3.51) 



d{Kpa'^) 
Jt 



-k{p + 3p)ad. 



aa 



1 d{Kpa'^) 1 



Aad = 0, 



6 dt 3 

el primer termino se puede escribir como la derivada del cuadrado de la derivada de a{t) 



1 d .2 1 d{Kpa'^) A d 



2 dt 6 dt 
escribiendo todo como una sola derivada 

d " 
'dt 



6 dt 



= 0, 



pa 



1 



0, 



3"^ 3 

por lo tanto, lo que esta dentro del parentesis debe ser igual a una constante —k 



• 2 ^ 2 

a — -Kpa 



1 



-Aa^ 



dividiendo por tenemos 



K 

3^+3 



A k 



(3.50) 



(3.51) 



(3.52) 



(3.53) 

(3.54) 
(3.55) 

(3.56) 

(3.57) 

(3.58) 
(3.59) 



en donde hemos usado de nuevo la definicion del parametro de Hubble H = ^. Vemos que podemos 
recobrar las ecuaciones dinamicas estandares de los modelos de FLRW a partir de la EDG. 



^Escribimos esta ecuacion multiplicando a ambos lados por k = SivG, la cual no quita generalidad a la expresion y se 
usa solo por estar de acuerdo con nuestras unidades. 
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3.3.2. EDG para Campos Vectoriales Nulos 



La EDG para campos vectoriales nulos nos dara de forma elegante las ecuaciones diferenciales para las 
distancias, pues considera la separacion de geodesicas nulas (rayos de luz o fotones). Debemos considerar 
la EDG para campos vectoriales nulos dirigidos al pasado. En este caso = k°', kak°' = 0, < 0. La 
ecuacion de Pirani (3.43) ahora es 



^(S-yS'^' 'I f^ — T.'^KP^F)^ '/ > (3.60) 

que se denomina Enfocamiento de Ricci. Asi, si escribimos r/" = rje", Cae" = 1, eau" = Cak" = y 
usando una base alineada y propagada paralelamente, ^ = A;^V^e" = 0, por la EDG encontramos 



1 



(3.61) 



En este caso, todas las geodesicas nulas dirigidas al pasado o al futuro experimentan enfocamiento si 
K.{p + p) > (mientras que el signo de A no tiene importancia). 



Se pretende obtener la ecuacion (3.61) escrita en terminos del parametro no affn redshift z. Veamos 
como se transforman los operadores diferenciales 

dz d 



d 



dv du dz ' 



(3.62) 



dz d 


(- 


dv dz 


ydu 




-2 - 


m 





\dz J dz"^ dz dz"^ 



(3.63) 



Para geodesicas nulas tenemos que 

(1 + ^) 



a 



E 

Eq 



dz 



1 + z 



da 
a 



(3.64) 

con a el factor de escala, y oq = 1 el valor presente (oq = 1. Entonces para el caso dirigido al pasado 

(3.65) 



dz 



,1 da ^ 
1 + z)-—di^ 
a dv 



[1 + z)-Edv = EoH(l + zY du, 
a 



donde de nuevo usamos la definicion para el parametro de Hubble. Asi obtenemos 

du 1 



dz EqH{1 + zY' 



dH 



expresamos como 



dz2 EqH{1 + z) 



dH dv dt dH 



\ , dH 



dH 



dz dz dv dt H{1 + z) dt ^ 



(3.66) 
(3.67) 
(3.68) 
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el signo menos viene de la condicion de geodesicas nulas dirigidas al pasado, cuando z aumenta, v 
disminuye. Usamos tambien que ^ = -EqII + -2)- Ahora, de la definicion del parametro de Hubble H 
tenemos 

if 



da d ^2 
dt a a ' 



de la ecuacion de Raychaudhuri (3.50) se obtiene 

d^i/ _ _ 3 

dz^ ~ ~EoH{l + z)^ 

y 

d'^rj 



1 + 



{p + 3p) 



+ 



l + z 

De esta manera, la ecuacion (3.61) se convierte en 

3 



1 + 



iEH{l + z)y 



d'^rj 
d^ 



l + z 



K A \ dr] 



9H^ / dz 



dividiendo por E'^ y organizando terminos 



d'^T] 3 
dz^ ^ l + z 



K , , A 

1+ TT7T777(p + 3p) 



+ 



1 



drj 

dz ' 2(1 + Z)2if2 



{p + p)v = 0. 



Si consideramos de nuevo una mezcla no interactuante de materia y de radiacion 

Kp = m^n^oii + zf + m^nroii + zf, kp = H^nroii + Z)\ 

usemos ademas la ecuacion de Friedmann 

= Hi [0^0(1 + zf + VLro{l + zf + + 17fc(l + zf\ , 
de modo que la EDG para campos vectoriales nulos es 

A 417^o(l + z)^ + (7/2)17^o(l + z)3 + SQkoil + zf + 217a dri 

dz"^ (1 + z) {Qroil + zY + 5^m0(l + zf + Ofco(l + zf + J7a) dz 

20,0(1 + ^)^ + (3/2)0^0(1 + ^) 



(3.69) 

(3.70) 
(3.71) 

(3.72) 

(3.73) 
(3.74) 
(3.75) 



+ 



r/ = 0. (3.76) 



Q,o(l + zY + rj^o(l + zf + J7fco(l + zf + 

1 — ^ImO — ^ro con lo cual 



La relacion de Mattig es obtenida en el caso Qa = y escribiendo Qko 
obtenemos [35] 



d'^r] 6 + a^oil + 7z) + ^'0(1 + 8z + \z^) drj 



dz^ 



+ 



2{l + z){l + Qmoz + nroz{2 + z)) dz 



2{l + z){l + nm0Z + nr0z{2 + z))'^ • I - J 



Vamos a considerar ahora la relacion entre la magnitud del vector desviacion rj con la distancia diametral 
angular. En un espacio-tiempo esfericamente simetrico, como el modelo FLRW, la magnitud del vector 
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desviacion rj se relaciona con el area propia dA de una fuente a redshift z por dr] oc V dA, y de esta 
expresion, la definicion de la distancia diametral angular d^ puede ser escrita como 




dA = \h?^, (3.78) 



con dil. el angulo solido. De modo que la forma de la relacion de Mattig es la misma cambiando ij por dA 
(el factor de proporcionalidad se puede descartar de la ecuacion diferencial) . Asi tendremos la siguiente 
ecuacion diferencial para la distancia diametral angular 

d'^dA 6 + nmo{l + 7z) + nro{l + 8z + 4z2) ddA 
dz'^ ~^ 2(1 + z){l + nmoz + ^rozi2 + z)) dz 

I 31^^o + 40,o(l + z) ^ 

2{l + z){l + nm0Z + nr0zi2 + z)) ^ ' 

Fijando entonces los parametros cosmologicos y resolviendo esta ecuacion diferencial, podemos encontrar 
todas las expresiones para las distancias cosmologicas, usando para esto las relaciones existentes entre 
ellas (3.30), (3.33). 
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Capitulo r 

Ecuaciones de Campo en Teorias f{R) 



The next case in simplicity includes those manifoldnesses in 
which the line-element may be expressed as the fourth-root 

of a quartic differential expression" 
Bernard Riemann, 1854 



unque la Relatividad General se supone valida, y se han probado varias soluciones experimental- 



mente, desde su genesis se han desarrollado teorias alternativas. En particular, Weyl y Eddington 
en 1919 propusieron teorias de gravedad que consistian en lagrangianos cuadraticos (proporcionales al 
cuadrado del escalar de Ricci), y se han estudiado las soluciones para estos espacio-tiempos, en parti- 
cular la de Schwarzschild [43]. 

El punto de partida fundamental de estas teorias de gravedad modificada consisten entonces en los 
lagrangianos y hamiltonianos que describen el campo gravitacional. La primera accion como vimos es 
la de Einstein-Hilbert 



como mostramos en el Capitulo 2 las ecuaciones de campo de Einstein se encuentran a partir de un 
principio variacional de esta accion. La accion para una teoria de gravedad cuadratica puede entonces 
escribirse como [44] 



en donde a es una constante con que permite tener el lagrangiano en la accion con las unidades correc- 
tas. Esta accion con un lagrangiano cuadratico conduce a unas ecuaciones de campo que son de cuarto 
orden en la metrica, lo que conlleva a que surjan nuevos grados de libertad en la teoria. Puede ademas 
pensarse en teon'as con lagrangianos compuestos con todos los posibles escalares que pueden obtenerse 
a partir del tensor de Riemann y sus derivadas, sin embargo tales teorias poseen el mismo problema de 
grados adicionales. 





(4.1) 
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Despues de que Einstein propusiera la RG y Hilbert encontrara el lagrangiano que la describe, Kretsch- 
mann [45] propuso una accion construida con el escalar R^^^y^R'^'^^^ en vez del escalar de Ricci, el cual 
fue denominado tiempo despues como el escalar de Kretschmann 

Sk= [ d^xV^R^r^^sR"''^^. (4.3) 
Jv 

Puesto que el tensor de Riemann Rari'yS es el tensor fundamental para la gravitacion y el escalar de 
Kretschamn es su cuadrado, este lagrangiano parece ser una buena opcion; ademas, esta teon'a preserva 
las identidades de Bianchi. Sin embargo, la aparicion de terminos de la forma V^V^i?"^^^, conduce a 
fijar los elementos de la metrica, sus primeras, segundas y terceras derivadas, de modo que la teorfa 
introduce mas grados de libertad en comparacion con RG. La aparicion de estos terminos de cuarto 
orden en la metrica se puede evitar considerando el escalar de Gauss-Bonnet [46] 

g = R^- ARapR""^ + RafSysR''^''^ (4.4) 
de modo que la accion pueda escribirse como 

Sg= [ d^x^g. (4.5) 
Jv 

Existe sin embargo una teorfa mas general de gravedades modificadas denominadas Teorias de Lovelock 
en la cual se da una generalizacion del escalar de Gauss-Bonnet [47], [48]. Tales teon'as se basan en un 
lagrangiano general de la forma 

r r — ^ri7,'S~^ JLo- ^1^^ P2--I32n 7j>/3l /32 . . . o/32n-l P2n (A (\) 

J-L — V y 2" " ')'lT2-'r2n -^7172 -^72n-l72n' y^'^! 

n 

en donde 5^l^2-'-'-^2n 1^ delta de Kronecker generalizada de orden 2n, y a„ son constantes. 

Otra posibilidad interesante, pero que despues de su aparicion fue en cierto modo olvidada, son las 
teorias siendo C la traza del tensor de Weyl, y denominadas teorias de gravedad conforme, que 
tienen la particularidad de ser una teorfa localmente invariante, y una de las candidatas para la formu- 
lacion cuantica de la gravedad. Las ecuaciones de campo en esta gravedad conforme fueron publicadas 
primeramente en [49] y su estudio se justifica en gran medida por ser una alternativa ffsicamente viable 
a problemas que la cosmologfa estandar posee, en particular la hipotesis de la materia oscura. 

Algunas extensiones de la accion de Einstein-Hilbert a modelos con escalares de orden superior son 
consideradas en [50]- [52]. Como una extension natural de RG y de las teorfas de orden superior surgen 
las teorias f{R) al considerar una funcion arbitraria del escalar de Ricci en vez de solamente el esca- 
lar como en el caso de la accion de Einstein-Hilbert [11]. En las teorfas de gravedad modificada f{R) 
debemos distinguir 3 casos diferentes, los cuales nos van a conducir a diferentes ecuaciones de campo. 
Est as 3 versiones se pueden resumir en 
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1. Formalismo Metrico: Para este caso escribimos la accion como 

Smod = Smet + S M {g aji ■, i^) ■, Smet = 1^ [ d'^X^/^f{R), (4.7) 

con K = SttG. Aqm la variacion se hace con respecto a g"'^^ y ^/^ denota todos los campos de 
materia. 

2. Formalismo de Palatini: En este caso la metrica conexion P^^ son independientes. Sea 
entonces TZa/s el tensor de Ricci construido con la conexion independiente y TZ = g'^^TZap- La 
accion puede escribirse como 

Spai = ^ I d^x^gfiTZ). (4.8) 

La accion total incluye la accion asociada a todos los campos de materia, la cual es independiente 
de la conexion. 

3. Formalismo Metrico-Afm: Esta forma es similar a la de Palatini, pero en este caso la accion de 
materia depende de la conexion 

Sma = ^ I d^x ^gfilZ) + SM{9ap,T%,^l,). (4.9) 

A continuacion mostraremos la obtencion de las ecuaciones de campo en los 3 formalismos, dando 
especial atencion a las ecuaciones de campo en el formalismo metrico, incluyendo terminos de frontera. 

4.1. Ecuaciones de Campo en el Formalismo Metrico 

Las ecuaciones de campo en el formalismo metrico fueron deducidas en [53], e incluyendo terminos de 
frontera en gravedades de cuarto orden en [54]. El termino similar al de Gibbons- York-Hawking en 
teorias f{R) fue explorado en [55], con un principio variacional aumentado en [56], y usando un marco 
tensor-escalar en [57]- [60]. 

Consideramos de nuevo el espacio-tiempo como un par [M^g) con Ai una variedad cuadri- dimensional 
y dap una metrica sobre M. El lagrangiano ahora es una funcion arbitraria del escalar de Ricci 
'^[g'a/3] = f{R)i de nuevo asumimos que la variedad posee una conexion de Levi-Civita, es decir, un 
simbolo de Cliristoffel. La accion general que incluye un termino de frontera se puede escribir como [60] 

Smod = -^{Smet + S'qy h) + ^M-, (4-10) 

Ik 

con 

Smet= I d^xV^f{R), (4.11) 

Jv 

y el termino del tipo Gibbons- York-Hawking [60] 

S'GYH = ^i d^yeyfi\f'{R)K, (4.12) 

JdV 
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con f'{R) = df{R)/dR. De nuevo, Sm representa la accion asociada con todos los campos de materia 
(2.42). Vamos entonces a encontrar las ecuaciones de campo en el formalismo metrico de gravedad f{R) 
siguiendo [61]. Realizamos la variacion con respecto al tensor metrico sujeto a la condicion 



0. 



dV 



Primero realizamos la variacion del termino Smet 

5Smet= [ d^x{f{R)5^g + ^g5f{R)), 
Jv 

la derivada funcional de f{R) se puede escribir como 

6fiR) = f'iR)6R. 
Usando la expresion para la variacion del escalar de Ricci 

6R = 6g"^Ra(s + V ^g"^ {ST'^p,) - ^"'^(^n^)), 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



en donde la variacion del termino g'^^{ST'^^) — g°'^{6rZ'y) esta dada en Apendice A. Con este resultado 
tendremos 

5R = dg'^^R^p + V45"''(<5r^J - g'^'i^Tl^)), 
= Sg'^^Rap + gf.uVaV'iSgn - V^V^iSg'^^), 

= dg'^^R^p + go.pUidg'^P) - VaVp{6g''^). (4.17) 

Con □ = Vo-V^ y renombramos algunos indices mudos. Reemplazando estos resultados en la variacion 
de la accion (4.14) 

SS„^et= [ d^x{f{R)5^g + ^gf{R)5R), 
Jv 

= ^ d^x (^-'^f{R)^g^^6g"P + f'{R)V^(^6g^f'R^^ + gafs^idg^^) - V,V^(55°^))) , 
= l^d^x^ (^fiR) (dg^^R^p + gapn{5g''^) - V^VpiSg''^)) - \f{R) g^pSg''^^ • (4.18) 
Consideramos ahora las siguientes integrates 

/ d''x^nR)g^pn{5g^P), [ d'x V^f iR)VaV p{6g"^). 
Jv Jv 



(4.19) 



Vamos a escribir estas integrates de una forma diferente realizando para esto una integracion por partes. 
Definimos las siguientes cantidades 



Mr = f{R)gapVr{Sg''^) - Sg'^^ go^pVrU'm, 
= nR)V^{5g'''<) - 5g^''V^if'{R)). 



(4.20) 
(4.21) 



4.1. Ecuaciones de Campo en el Formalismo Metrico 
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La combinacion g'^'^Mr + A^*^ es 

g'^^Mr + N^ = nR)g^pV%6g^^) - 5g^^ g^pV%f {R)) + f {R)V ^{6g"^) - 5g"^V.,{f {R)), (4.22) 

en el caso particular f{R) = R, la anterior combinacion se reduce a la expresion (2.24) con la ecuacion 
(A. 10). Las cantidades Mr y nos permiten escribir la variacion del termino (4.18) de la siguiente 
forma (ver Apendice B) 

+ 1 cPy£y/\h\n'Mr+ I d^y ey/\h\n„N'' . (4.23) 
JdV JdV 

A continuacion estudiaremos la contribucion en la frontera (4.23), y mostraremos que este termino se 

cancela con las variaciones de la accion Si- 



GYH- 



4.1.1. Termino de Frontera en el Formalismo Metrico f{R) 

Expresamos las cantidades y calculadas en la frontera dV. Es conveniente primero expresar estas 
cantidades en funcion de las variaciones Sgap- Usando la ecuacion (2.15) en (4.20) y (4.21) obtenemos 

Mr = -nR)g"Nr{6g^p) + g^^^gap^ rU'iR)). (4.24) 



= -nR)g''^g<''y-,{5g^u) + g'^'^f'^dg^^V^ifiR)). (4.25) 

Para evaluar estas cantidades en la frontera usamos el hecho que 6gai3\dv = Sg^^ldv = 0; de modo que 
los linicos terminos no nulos son las derivadas de 5gai3 en las derivadas covariantes. Asi tendremos 



Mr 



dV 



dV 



-fiR)g''^dr{dg^p), 
-f'{R)g^^g'y'd^{6g,,), 



(4.26) 



(4.27) 



Calculamos ahora n'^Mrlgy y no-A^'^j^^ que son los terminos en las integrales sobre la frontera (4.23) 



lav 



n'^Mr 



dV 



-nR)n-'{en^n^ + h^P)dr{5g^p), 



= -f'{R)n''h''^d^{5ga 
donde renombramos el indice mudo r. Por otro lado 



(4.28) 



dV 



-f'iR)n^{h'^'' + en'^nnd^i6g^u), 
-f{R)n^'h^''d^{6g^,), 



(4.29) 
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donde usamos que Ua h''^' = 0, = 1 y el hecho de que la derivada tangencial h^'^ d-y{6g^y) se anula. 
Con estos resultados la variacion de la accion Smet es 

- / d'ye^/lhlf'iRKh'^f'd^iSg^p). (4.30) 
JdV 

Procedemos ahora con el termino de frontera Sqyh accion total. La variacion de este termino es 

SS'cYH = 2 / e^\{6f{R)K + f'iR)6K) , 

= 2 I d^yey^\{f"{R)6RK + f'{R)6K). (4.31) 

JdV 

Usando la expresion para la variacion de K, ecuacion (2.40), podemos escribir 

SS'gyh = 2 £^ d^ye^\(^f"{R)5RK+ \f {R)h^^d^{5gp^)n''^ , 

= 2 I d^y£y^\f"{R)6RK+ (f d%ey^\f'{R)h''Cd^{6gpa)n''. (4.32) 
JdV JdV 

Observamos que el segundo termino en (4.32) cancela el termino de frontera en la variacion (4.30), y 
adicionalmente necesitamos imponer 5R = en la frontera. Un argumento similar esta dado en [60]. 

Finalmente, con la variacion de la accion de materia, dada en (2.45), la variacion total de la accion en 
el formalismo metrico de gravedad f{R) es 



5S, 



mod 



-^j^d''x^To,p6g''^. (4.33) 



Imponiendo que esta variacion sea estacionaria tendremos 

1 SSmod 



-9 Sg^P 



=^ f{R)RaP + gafsOfiR) - VaV(sf'{R) - -f{R) gafS = 



de modo que las ecuaciones de campo en el formalismo metrico f{R) son 



f'{R)Ra^ - ^f{R) ga/3 + QapUfiR) - V,V^/'(i?) = KT^p. 



(4.34) 



(4.35) 



Estas ecuaciones son claramente ecuaciones de cuarto orden en la metrica (por la contribucion de 
segundas derivadas en R de los operadores □ y VqV/j). La traza de estas ecuaciones de campo es 



f{R)R - 2f{R) + SnfiR) = kT, 



(4.36) 



en donde podemos observar que a diferencia de la RG, la relacion entre R y T pasa de ser una simple 
relacion algebraica a una relacion diferencial para R. Este hecho conlleva a resultados interesantes. 



4.2. Ecuaciones de Campo en Formalismo de Palatini 
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El teorema de Birkhoff que asegura que la solucion de Schwarzschild es la solucion mas general a las 
ecuaciones de campo en el vacio, deja de ser valido en teorfas f{R)- En RG la relacion entre R y T 
viene dada por R = —kT, lo que implica que en el vacio las ecuaciones se reducen a, R = 0; en este caso 
sin embargo si T = no necesariamente se obtiene un escalar nulo. 



Es posible escribir las ecuaciones de campo en la forma de las ecuaciones de Einstein con un tensor 
de energia-momentum efectivo compuesto de terminos de curvatura. Para esto escribamos la expresion 
(4.35) en la forma 

1 



R 



iRgalz 



kT^P ^ [f{R) - Rf'jR)] ^ [VaV^fiR) - gai^afiR)] 



f'{R) 



2f>{R) 



f'{R) 



G. 



f'{R) 



donde 



rpeff _ 



[f{R)-Rf'm 



9al3 + [VaVp-gaf}a]fiR) 



(4.37) 
(4.38) 

(4.39) 



De esta manera los efectos de la modificacion en las ecuaciones de campo se pueden interpretar como 
una contribucion de un fluido de curvatura correspondiente a terminos netamente geometricos. Se pue- 
de mostrar que este tensor de energia-momentum efectivo satisface tambien una ley de conservacion 



A continuacion mostraremos muy resumidamente las ecuaciones de campo en los formalismos de Palatini 
y metrico-afm siguiendo [63] , [64] . 



4.2. Ecuaciones de Campo en Formalismo de Palatini 

La idea basica en este formalismo es que el tensor metrico gajs y la conexion P^^ son independientes. 
Consideremos entonces TZa/s el tensor de Ricci construido a partir de esta conexion independiente y sea 
TZ = g'^^TZaji el respectivo escalar de Ricci. Partiendo de la accion dada en la ecuacion (4.8), y haciendo 
las variaciones con respecto a g"^^ y P^^, respectivamente, encontramos las ecuaciones 

S'{n)n(^o.p) - \f{n) g^p = kT.^j, (4.4o) 

- V^(^/=5/'(7e)5"^) +V.(V^/'(7^)5'^(")5^) = 0, (4.41) 

con la derivada covariante definida con la conexion independiente, y en donde los indices entre 
parentesis estan simetrizados. Tomando la traza de la ecuacion anterior se pueden reducir las ecuaciones 
a 

/'(7^)7^(,^)-^/(7^)g«^ = KT,;3, (4.42) 
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V,{^gnn)g^^)=^. (4.43) 
Consideremos la ecuacion (4.43), definamos una metrica conforme a g^p de la siguiente forma 

Kp = f{n)g^p, {AAA) 

por lo tanto la ecuacion (4.43) se escribe como 

V^(V^/i"'') = 0, (4.45) 

que se puede interpretar como la conexion de Levi-Civita para /iq,^ 

T% = h^'^idaKf^ + d^Ka - d^Kp) (4.46) 
= ^-^5"n^«(/'(^)5./3) + dp{f{n)g,^) - dM'(.T^)9a(s)). (4.47) 

Puesto que la traza relaciona TZ algebraicamente con T, y como tenemos una expresion explicita para 
en terminos de 7^ y 5"^, podemos en principio eliminar la conexion independiente de las ecuaciones 
de campo y expresar esta solo en funcion de la metrica y los campos de materia. Veamos ahora como 
se escribe el tensor de Ricci bajo una transformacion conforme 

7^a/3 = Rap + lj^;^{^M'm)i'^pifm) - 77^(V"V/3 - 9alSa)f\n), (4.48) 

contrayendo esta expresion con g"^ nos da 

^ = ^ + ^^^^;^(v.(/'(^)))(v"(/'(7^))) + (4.49) 

con lo cual podemos obtener 



Gaf^ = j^fa, -l{n- j^)9a, + (v.v, - lg^,n^ fin) 



^ ^ ^{v^f'{n)){v^f'{n)) - \g^p{vf\n)f ] , (4.50) 



2(/'(7^))2V' ' ' " " 2 

en donde escribimos 

Gap = Rap — -j^Rgap- (4-51) 

Como es de esperar, las ecuaciones de campos en los formalismos metrico y de Palatini se reducen a las 
ecuaciones de campo de Einstein en el caso f{R) = R, o con la contribucion de la constante cosmologica 
cuando f{R) = R-2A. 



4.3. Ecuaciones de Campo en el Formalismo Metrico-Afm 
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4.3. Ecuaciones de Campo en el Formalismo Metrico-Afm 

Este formalismo es similar al de Palatini, excepto que ahora los campos de materia dependen de la 
conexion F^^. Ademas se esto se da la libertad para que la variedad A4 en este formalismo tenga 
torsion (F^^ 7^ -'^7/3 )• Esta libertad permite que en algunos regi'menes de energia el espm de las parti'culas 
pueda interactuar con la geometria, lo cual directamente nos da la existencia de la torsion en la teoria. 
Considerando entonces la accion (4.9), y usando los siguiente resultados para la variacion del escalar TZ 
(asumiendo torsion diferente de cero) 



(4.52) 



y la variacion de la accion de materia Sm = 'S'a/(5q:/3 5 T^^, V') 



55m ^ 7 



donde definimos ahora un nuevo tensor llamado el Hipermomentum como 



7 



2 5Sm 

al3 



n 



llegamos a las ecuaciones de campo 



(4.53) 



(4.54) 



/'(7^)7^(a/3) - -fi^R) gap = KTa 



(4.55) 



+ 27=g /'(7^) {g'^^S!^^ - g^PS;, + g'^'^S^^^) = ^Af , (4.56) 

donde escribimos la ultima expresion en funcion del tensor de torsion de Cartan (2.3). Un ejemplo de 
campos que tengan hipermomentum diferente de cero, son los campos de Dirac, asi que el estudio de 
tales campos en este formalismo nos podrian dar evidencia de la existencia de torsion en la teoria [64]. 



4.4. Equivalencia con Teorias de Brans-Dicke 

Una formulacion diferente a la RG que ademas del tensor metrico (y de este los tensores de curvatura) 
incluye un campo escalar ip, fueron trabajadas por Brans y Dicke en 1961 [65]. La motivacion principal 
de estas teorias era la de cambiar el valor de la constante de gravitacion universal G de lugar a lugar, 
con el fin de incorporar el principio de Mach^ en la RG. 

^Que en resumen nos dice que la inercia de un cuerpo depende en principio de todo el contenido de materia del universe. 
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La accion general para las teon'as de Brans-Dicke se puede escribir como 



S 



d X \/— (7 



{Vct>f - U{ct>) 



+ SM{ga(5,tp), 



(4.57) 



en donde lobd es el denominado parametro de Brans-Dicke y (V(/))^ = g^'^da'pd/scl)- 



Las teorias f{R) en el formalismo metrico se pueden escribir en la forma de las teorfas de Brans-Dicke 
[11] con un potencial para el campo escalar como grado de libertad (escalaron). Para esto consideramos 
la siguiente accion (sin terminos en la frontera) para un nuevo campo x 



2k 



(4.58) 



en donde /^^ corresponde a la derivada de la funcion / con respecto a cabe aclarar que el termino 
f^xi^ ~ x) 6S un producto entre f^^ Y x)- Variando esta accion con respecto a x obtenemos 



lxxiR-x) = o, 



(4.59) 



fijando f^^x 7^ se sigue que i? = x- De este modo la accion (4.58) recupera la accion (4.7) en teorfas 
f{R)- Si definimos ahora 

<t>^f,xix)^ (4-60) 

la accion (4.58) puede escribirse como 



S 



donde U [(f) esta definido como 



d X ■ 



2k 



(pR - U{(p) 



+ SA4{gai3,i'), 



U{ct>) 



x{4>)^- f{x{(t>)) 



2k 



(4.61) 



(4.62) 



De esta forma si comparamos la expresion (4.57) con (4.61) observamos que las teorias f{R) con equi- 
valentes a las de Brans-Dicke con wbd = (y en unidades de k = 1). En el formalismo de Palatini 
donde la metrica gaj3 y la conexion F^^ con considerados como variables independientes, el escalar 
de Ricci es diferente que el dado en el formalismo metrico. Se puede mostrar que para el formalismo 
de Palatini en teorfas f{R) es equivalente a las teorfas de Brans-Dicke con el parametro wbd = —3/2 [11]. 



Hemos entonces discutido los principios variacionales en teorfas f{R), haciendo especial enfasis en 
el formalismo metrico. Nuestro punto de partida entonces para construir modelos cosmologicos son las 
ecuaciones de campo modificadas (4.35). 
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Capitulo 
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Modelos Cosmologicos en Teorias f{R) 



-rj\ n este capitulo aplicaremos las ecuaciones de campo en el formalismo metrico de las teorias f{R) pa- 
ra un universo homogeneo e isotropico descrito por la metrica de Robertson- Walker. Varias formas 
de la funcion f{R) se han propuesto con el fin de obtener la actual expansion acelerada del universo, 
sin embargo, algunos poseen problemas (en particular inestabilidades) en el regimen de bajas energias 
(sistema solar). 

El estudio de de las ecuaciones de Friedmann modificadas, modeladas como un sistema dinamico au- 
tonomo, permiten restringir la forma de las funciones y dar criterios de viabilidad para teorias que 
reproduzcan una epoca de radiacion, una de materia y una de expansion acelerada. 



5.1. Ecuaciones de Friedmann modificadas en gravedad f{R) 

Debemos entonces hallar las expresiones para las ecuaciones de Friedmann modificadas asumiendo para 
esto de nuevo la validez del Modelo Estandar de la Cosmologfa. De modo que nuestro punto de partida 
es la metrica de Robertson- Walker, ecuacion (3.1) 



-dr + a\t) 



+ r'^de'^ + sin^ 



1 — kr"^ 



(5.1) 



siendo a{t) el factor de escala y k \a. curvatura espacial del universo. Para esta metrica el escalar de 
Ricci es 



i? = 6 



a f d\'^ k 


= 6 




a \a J 





H + 2H^ + 



k 



(5.2) 



consideraremos tambien el tensor de energia-momentum de un fluido perfecto, ecuacion (3.2) 

To,fi = {P + P)UaUl3 + pgal3, (5.3) 

siendo p la presion del fluido, p la densidad de energfa y u" la cuadrivelocidad de los observadores 
fundamentales. A partir de las ecuaciones de campo (4.35), y con la metrica (3.1), obtenemos las 
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siguientes ecuacion de Friedmann modificadas (ver Apendice C) 



3f'{R) 



2H + 3i?^ + 



Kp + 2HRf"{R) + (/(-^) Rf'iR)) ^ ^^//(^) ^ {Rff"{R) 



(5.4) 



(5.5) 



en donde usamos que H = ^. A partir de la traza de la ecuacion de campo (4.36), podemos encontrar 
otra relacion para f{R) y sus derivadas con la presion y la densidad de energia 



f{R)R - 2f{R) + 2,Uf{R) = 43p - p). 



(5.6) 



A partir de las ecuaciones (5.4) y (5.5), podemos definir una densidad y presion asociadas a los terminos 
adicionales de curvatura de la siguiente forma 



1 



Pde 



{Rf'iR)-fiR)) 



3Hf"{R)R 



PDE 



2HRf"{R) + 



if{R) - Rf'm 



+ Rf"{R) + {RYf"{R) 



de modo que la ecuacion de estado puede escribir se como 



WDE 



PDE _ 2HRf"{R) + (f(R)-^f'(R)) + Rf{R) + {Rff"'{R) 



PDE 



-1 + 



(Rf'(R)-f(R)) _ 2,Hf"{R)R 

Rf"{R)-HRnR) + {R)'r{R) 
(Rf'{R)-f{R)) 



2 3Hf"{R)R 

De forma similar a la ecuacion de continuidad para las contribuciones de densidad y materia 

Pm + 3Hpm = 0, 



(5.7) 
(5.8) 



(5.9) 



(5.10) 



Pr + 4:Hpm = 0, 



(5.11) 



podemos obtener una expresion para la ecuacion de estado efectiva para las componentes poE Y Pde 
(ecuaciones (5.7) y (5.8) respectivamente) [7] 



PDE + 3-^(1 + Wde)PDE = S-ffo^m-^^y-^^ 

(/ \^>) 



(5.12) 



Si consideramos ahora las ecuaciones de campo modificadas constituidas por un fluido efectivo (4.39), 
la ecuacion de Friedmann (3.4) se puede escribir como 



a K 
a 6 



Ptot + 3ptot), 



(5.13) 
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en donde el submdice tot denota la suma del fluido de curvatura y la contribucion de materia a la 
densidad de energia y presion. De esta relacion la condicion de aceleracion, para un modelo en el 
dominio de materia es 

PDE + Pm + 3pDE <0, WDE<--^^, (5.14) 

PDE 

esta relacion nos permite entonces determinar si dada alguna forma especifica de la funcion f{R) y por 
lo tanto Unas componentes para poE y Pde, obtener universos con expansion acelerada. A los modelos 
que tienen una ecuacion de estado con w < —1/3 los denominaremos modelos phantom o de Quintae- 
sencia. 



Aunque algunos modelos se han utilizado para explicar la expansion del universo, en particular el 
modelo [10], [66] 

ct 

f{R) = R-—, (a>0,n>0), (5.15) 

presenta sin embargo un numero de problemas, en particular la inestabilidad en la materia [18], [67], 
inestabilidad en perturbaciones cosmologicas [68]- [71], y la incapacidad de satisfacer restricciones locales 
(sistema solar) [72]- [74]. Este tipo de pruebas han llevado a considerar algunas restricciones en las 
funciones f{R) [21] en particular 



f{R) > 0, f"{R) > para R > i?o(> 0), (5.16) 



siendo Ro el valor actual para el escalar de Ricci. Consideraremos a continuacion la aproximacion 
de sistema dinamico en teorias f{R) que nos permitira obtener mas criterios de viabilidad para las 
funciones f{R). 



5.2. Aproximacion de Sistema Dinamico en ecuaciones de Friedmann 
modificadas 

El analisis de las ecuaciones de campos modificadas modeladas como un sistema dinamico, nos permitira 
encontrar criterios para modelos cosmologicamente viables. Realizaremos un analisis de tal enfoque 
siguiendo [ll],[75]-[77]. Vamos primero a escribir las ecuaciones de Friedmann modificadas (5.4) y (5.5) 
en el caso de un universo piano (fc = 0) y considerando las contribuciones en la densidad de energia y 
la presion para los dominios de radiacion y materia. Para esto usamos que la ecuacion de estado para 
el caso de materia es pm = y en el caso de radiacion pr = l/3pr, con esto las ecuaciones pueden 
escribirse como 

3f'H^ = K{pm + Pr) + - 3Hf', (5.17) 

-2f'H = Jp^ + ^p,\+f'-Hf\ (5.18) 



5.2. Aproximacion de Sistema Dinamico en ecuaciones de Friedmann modificadas 
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en donde usamos la siguiente notacion /' = dt{f') = f"R, y eliminamos la dependencia de R en las 
funciones / y /'• Con M = lna, definimos un operador diferencial de la siguiente manera 



d d Id 



dM dlna H dt' 
asi que la ecuacion (5.20) se puede escribir como 



(5.19) 



^ Kpm Kpr _ / _ , . 

Definamos las siguientes cantidades adimensionales 

XI = -^, ^2 = -g^, X3 = 2 + -^, X4 = 3^, (5.21) 

= ^jr^ = 1 - - X2 - X3 - X4, (5.22) 

y tambien 



^DE = Xi + X2 + X'i, (5.23) 



el cual lo asociaremos con la contribucion adicional a las ecuaciones de campo, o equivalentemente la 
contribucion a la energfa oscura (razon por la cual usamos los submdices DE: Dark Energy). Con estas 
definiciones se pueden obtener las siguientes ecuaciones de movimiento 

= -1 - X3 - 3x2 + X1X3 + X4, (5.24a) 

j^ = — -X2(2x3-4-xi), (5.24b) 
dJ\ m 

dx-i X1X3 , , , , 

_5 = _^_2.3te-2), (5.24c) 

= -2x3X4 + X1X4, (5.24d) 

donde hemos definido 

_ din/' _ Rf" _ Rf _ din/ _ Rf _ X3 

"^-dlnR- f ~ f ' d\nR~ f " X2 ' ^ 

De la expresion para r, el escalar R se puede escribir como X3/X2, y puesto que m depende de R, esto 
significa que m es una funcion de r, m = m{r). El modelo ACDM, f{R) = R — 2A, corresponde a 
m = 0. Por lo tanto la cantidad m caracteriza la desviacion en la dindmica de un modelo ACDM. A 
partir de las ecuaciones (5.20) y (5.18) podemos definir una ecuacion de estado efectiva por la siguiente 
relacion 

2H (2X3 - 1) 
-e// = -l-^ = (5.26) 

la cual se puede interpretar como la ecuacion total para una contribucion de densidad de energfa efectiva 
Peff = Pm + Pr + Pde y Una relacion similar para la presion efectiva. 
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P : (xi,X2,X3) 



eff 



Pi = (0,-1,2) 



P2 = (-1,0,0) 



^3 = (1,0,0) 



P4 = (-4,5,0) 



ft 



3m- 



1 + 4m 1 + 4m 



1+m' 2(1 + m)2' 2(1 + m) 



2(1 -m) 1 - 4m -(l-4m)(l + m) 
1 + 2m ' m(l + 2m) ' m(l + 2m) 



1/3 



1/3 



1/3 



m 
1+m 



2 — 5m — 6m 
3m(l + 2m) 



m(7 + 10m) 
2(1 + m)2 



Tabla 5.1: Valores para los puntos cn'ticos del sistema dinamico (5.24) y valores para la ecuacion de estado efectiva We// 
(5.26). 



La idea ahora es determinar la estabilidad del sistema dinamico. Vamos a realizar el siguiente ana- 
lisis en el caso en que no hay radiacion, es decir, para X4 = 0. Notamos primero que las ecuaciones 
dinamicas se pueden escribir en la forma 

dx ' 

= gi{xi,X2,X3), i = l,2,3, (5.27) 

de modo que los puntos fijos del sistema son los puntos ( satisfacen de forma simultanea 

gi{x'i, X2, x'^) = gj{x'i,X2,x'^) = = 1,2,3, i 7^ j. En la tabla 5.1 se muestran tales puntos, el 

respectivo valor para ecuacion de estado efectiva (5.26)y el valor para el parametro Clm- Los puntos P5 
y Pq satisfacen la ecuacion 

X3 = -{m{r) + l)x2, (5.28) 

es decir 



m{r) = -r - 1. (5.29) 

Cuando m{R) no es constante, se debe resolver tal ecuacion, para cada raiz uno obtiene un punto 
del tipo P5 o Pg con m = m(rj). 

5.2.1. Estabilidad de los puntos criticos 

Para determinar la estabilidad de los puntos criticos vamos a seguir [78]. Consideremos por el momento 
el siguiente sistema de ecuaciones diferenciales acopladas para dos variables x{t) y y{t) 



x = f{x,y,t), y = g{x,y,t), 



(5.30) 
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donde f y g son funciones de x, y y t. Como vimos anteriormente, un punto cn'tico {xc,yc) es un punto 
tal que 



if ^9) 



0. 



(5.31) 



Para estudiar la estabilidad vamos a considerar pequenas perturbaciones 6x y 5y alrededor de los puntos 
cnticos {xc,yc), es decir 

X = Xc + Sx, y = y^ + 5y. (5.32) 

de modo que reemplazando esto en el sistema dinamico obtcnemos que 



d 5^ 



J 



(5.33) 



donde nuevo = In a y con J' dado por 



J 







dx 


dy 


dl 


dl 


\dx 


dy) 



(5.34) 



Esta matriz posee dos autovalores j\ y j2, de modo que la solucion general a la evolucion de las 
perturbaciones lineales es 

Sx = Cie^i^ + C2e^^^, (5.35a) 

Sy = Cge^i-^ + C4e^2-^, (5.35b) 

siendo Ci,C2,Cs,C4 constantes de integracion. Entonces, la estabilidad de los puntos fijos depende de 
la naturaleza de los autovalores 

■ (i) Punto Estable: ji < y j2 < 0. 

■ (ii) Punto Inestable: ji > y j2 > 0. 

■ (in) Punto de silla: ji < y j2 > (o ji > y j2 < 0) . 

■ (iv) Espiral estable: El determinante de la matriz J es negativo y las partes reales de ji y j2 son 
negativas. 

Un punto fijo es un atractor en los casos (i) y (iv), pero no lo es en los casos (ii) y (Hi). 



De esta manera construimos entonces la matriz J a partir de las derivadas parciales de las funcio- 
nes gi con respecto a cada una de las variables Xi. 

^2x1 — X3 —3 —1 — xi \ 



J(xi,X2,X3) 



— +X2 XlX3dx2\ — 

m \m 



\ m 



-XiX2,dx 



2x3 + 4 + Xl 



xidx-. 



2X2 



-xidx,{^ 1-4x3 + 4 



(5.36) 
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determinamos los valores propios del jacobiano a partir de la relacion 

det(^L,,,,,,3,-AI) = 0, (5.37) 

con I la matriz identidad. La estabilidad de los puntos criticos dependera entonces de los signos de los 
autovalores de la matriz (5.36) evaluados en los puntos criticos. 

• Pi: Punto de de Sitter. 

Puesto que Wgjf = — 1 el punto Pi corresponde a las soluciones de de Sitter que satisfacen 

F{R)R-2f{R) = 0, (5.38) 

que corresponde a la traza de las ecuaciones de campo (4.36) en el vacfo (T = 0). Para este punto los 
autovalores de (5.36) son 

donde mi = m{ri),ri = —2. Asi', Pi es estable cuando < mi < 1 y un punto de silla en los demas 
valores. La condicion entonces para la estabilidad del punto de de Sitter es 

< m{r = -2) < 1. (5.40) 



• P2: (j)-EDU. 

Este punto es caracterizado por una epoca "cinetica" en la cual la materia y un campo (p coexisten 
con fracciones de energia constantes. Denotamos entonces este punto por (/)-EDM (0-Epoca de Domino 
de Materia) siguiendo [79]. Los autovalores son 



7+- 



1712 rnl 



r(l+r)±W 7 + 



1 



nir, 



1712 fno 



r(l + r) -4 12 



m2 my 



r(3 + 4r) 



, (5.41) 



en donde m2 = (r2). Si m(r) es constante, los autovalores se reducen a 



4 + 



m2 



de modo que P2 es un punto de silla. 



(5.42) 



P3: Punto puramente cinetico. 



Este punto tambien corresponde a una epoca cinetica, pero a diferencia del punto P2 la fraccion de 
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energi'a en la materia se hace cero. Este punto ademas se puede ver como un caso particular del punto 
Pe escogiendo m = 1/4. Los autovalores son 



9+ — + ^r(l + r)±Wr9+ — -^r(l + r)V-4('20- — + ^r(5 + 4r) 



(5.43) 



Si m(r) es constante los autovalores se reducen a 



5, 2, 4 + —, (5.44) 

ma 

en este caso P3 es inestable para nis < and 1713 > 1/4, y un punto de silla para los demas casos. 

• P4. 

Este punto tiene una propiedad similar al punto P3 pues ambos tienen los mismos valores para w^ff y 
^rn- Tambien es un caso especial del punto Pq con m = —1. Los autovalores son 

4 

-5, -3 4+ — , (5.45) 
m4 

de modo que es estable para — 1 < m4 < y un punto de silla para los demas valores. Ninguno de los 
puntos P3 o P4 puede ser usado para la epoca de dominio de materia, ni para la epoca de expansion 
acelerada. 

• P5. 

El punto P5 correspondc a scaling solutions las cuales dan una razon constante Q„JO,£)e- En cl If- 
mite ms — > el punto rcprcsenta una era de materia estandar con a cx t^/'^ y = 1. De modo que 
la condicion necesaria para que P5 sea una era exacta de materia esta dada por 

m{r = -1) = 0. (5.46) 

Los autovalores para P5 son 

/. -3m5 ± A/m5(256m| + 160m| - 31m5 - 16) 

4m5(l + ms) 

en el li'mite 1 777,5! ^ ^ autovalores se pueden escribir aproximadamente por 



3(l + m'5), (5-48) 

para modelos con < 0, las soluciones no pueden permanecer por un periodo largo de tiempo alrededor 
del punto P5 debido al comportamiento divergente de los autovalores cuando m — > 0~ [m tendiendo a 
cero por la izquierda). El radical de los dos liltimos autovalores sc hacc negative para < < 0.327, 
de modo que los autovalores son complejos con partes reales negativas, asi que exigiendo que m'g > — 1, 
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(A) 




(B) 

A 


(D) 


(C) 



-3-2-10 1 2 3 



m, 

Figura 5.1: Comportamiento de We// en funcion de m para el punto Pe, obtenido de [76]. 

el punto P5 corresponde a un punto de silla. La condicion entonces para que este sea un punto de silla 
as 

m(r < -1) > 0, m'{r < -1) > -1. (5.49) 
• Pq: Punto de dominio de curvatura. 

Este corresponde a un punto de dominio de curvatura [76] cuya ecuacion de estado efectiva depen- 
de de m. La condicion para un universo con expansion acelerada {wefj < —1/3) se obtiene cuando 
mg < —(1 + \/3)/2, —1/2 < mg < 0, y itiq > (\/3 — l)/2. En la Figura 5.1 se muestra el comportamien- 
to de Weff en funcion de m Los autovalores para este punto son 

l-4m6 2-3me- 8m^ 2(1 - m6)(l + m^g) 
ruQ ' m6(l + 2m6) ' mQ{l + 2me) 

de modo que la estabilidad de Pq depende tanto de me como de m'g. En el Hmite mg — > ±00 tenemos 
Pq —7- (—1,0,2) con una ecuacion de estado w^jf = —1. Este punto es estable con m'g > —1. Pq es 
tambien un punto de de Sitter para uiq = 1, que coincide con Pi, y puesto que en este caso r = — 2, 
este punto se caracteriza por 

m(r = -2) — >1. (5.51) 
El punto Pq es estable y acelerado para 4 diferentes rangos: 

(I.) m', > -1 

Cuando uiq > —1, Pq es estable y acelerado en las tres siguientes regiones: 

(A) ruQ < —(1 -|- \/3)/2: Pq es acelerado pero no con condicion de phantom cosmology, es decir. 
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P : 


(xi,X2,X3) 








Weff 




Pi 


= (0,0,0,1) 








1/3 







/ 4m- 


2m 


2m 


1 — 2m — 5m^ \ 


1 — 3m 





\ 1 + m ' 


(l + m)2' 


(1 + m) 


' (l + m)2 ) 


3(1 +m) 



Tabla 5.2: Valores adicionales para los puntos cn'ticos del sistema dinamico (5.24) considerando contribuciones en radia- 
cion. 

Weff > —1. Se tiene el Hmite Weff = — 1 cuando uiq — )■ — oo. 

(B) — 1/2 < mg < 0: Pq obedece una ecuacion de estado phantom con Weff < —7.6. 

(C) rriQ > 1: Pq tiene una ecuacion de estado phantom con —1.07 < w^fj < —1. Se tiene 
Weff = —1 en el limite mg — ?• +oo y mg — )• 1. 

(II.) m', < -1 

Cuando mg < —1, Pq es estable y acelerado en la siguiente region: 

(D) {V3 + l)/2 < mg < 1: en este caso Pq no tiene una ecuacion de estado phantom Weff > — 1. 

De esto se deriva entonces una primer conclusion general con respecto a las funciones f{R) [7(5]: La 
aceleracion asintdtica NO puede tener una ecuacion de estado en el rango —7.6 < w^ff < —1.07. 

Si consideramos ahora el caso con radiacion los puntos Pi_g se mantienen iguales (con X4 = 0) y 
ademas obtenemos dos nuevos puntos, Tabla 5.2. Observamos que el punto Pj es un punto estandar de 
radiacion. Cuando m(r) es constante los autovalores de Pj son 

1, 4, 4, -1, (5.52) 

lo que implica que P7 es un punto de silla. El punto Pg es un punto nuevo de radiacion que contiene 
contribucion a la energi'a oscura diferente de cero. Puesto que la ecuacion de estado efectiva esta res- 
tringida por nucleosmteis a ser cercana a 1/3, es un punto aceptable para la epoca de radiacion con 
ms cercano a 0. Los autovalores para el punto Ps son 

En el hmite mg — > los liltimos dos autovalores son complejos con partes reales negativas, lo que 
muestra que el punto Pg es un punto de silla alrededor del punto de radiacion. A diferencia del punto 
de dominio de materia P5 no existen singularidades en los autovalores para el punto Pg en el hmite 
mg — > 0. Tambien notamos que el punto Pg esta sobre la hnea m = —r — 1 como en el caso del punto P5 
. Si la condicion para la existencia del punto de dominio de materia P5 se satisface (m(r) ~ 0~^,r = — 1), 
entonces existe un punto de radiacion en la misma region [76]. 

De este modo existen dos clases de modelos que son cosmologicamente viables: 



54 



Capi'tulo 5. Modelos Cosmologicos en Teon'as f{R) 



■ (A) Modelos que conectan -P5(r ~ —1, m ~ 0^) con Pi{r = —2, < m < 1)) 

■ (B) Modelos que conectan i-*5(r ~ — l,m ~ 0^) con PQ{m = —r — 1,— (\/3 + l)/2 < m < 0). 

De la expresion (5.16) los modelos viables de f{R) para explicar energia oscura deben satisfacer m > 
que esta en acuerdo con los argumentos anteriores. 



5.3. Algunos Modelos Especificos 

En esta seccion consideraremos algunos modelos f{R) para los cuales se puede escribir explicitamente 
m como funcion de r y se estudiara la posibilidad de obtener una epoca de dominio de materia seguida 
por una de expansion acelerada, todo esto siguiendo [76] 

5.3.1. f{R) = aR-"" 

Para este modelo obtenemos a partir de (5.25) la siguiente expresion 

m = —n — 1, (5.54) 

con r = n. La curva m(r) es degenerada, de modo que este caso se reduce a un sistema bidimensional 
(en el caso en que no hay radiacion), debido a la relacion xs = rx2 = nx2- La condicion m{r = —1) = 
se satisface entonces para n = —1, es decir, en el caso de la Relatividad General. Aliora para que el 
termino f'{R) = —naR^^^^ sea mayor que cero, se necesita que a < para n > y a > para 
n < 0. De la expresion (5.54) se tiene m > para n < —1 y m < para n > —1. Ahora a partir 
de los autovalores para el punto P5 (5.47), el punto para la epoca de materia es una espiral estable 
para n < —1, de modo que las soluciones no abandonan la era de materia para entrar en la expansion 
acelerada. 

Por otro lado, P5 es un punto de silla para — 1 < n < —0.713 mientras que el punto P2 es estable 
en la region —1 < n < —3/4. Sin embargo, uno de los autovalores de P5 muestra un divergencia positiva 
en el Hmite m — )■ 0~, lo que significa que el punto de materia se vuelve repulsivo para m cercano a 0~. 
Como vimos anteriormente la ecuacion de estado efectiva para Pg tiene un valor phantom de < —7.6 
para m = 0~, lo que implica que el punto de Silla P5 se conecta al punto P2 o al punto Pq. 

El modelo ACDM f{R) = R — 2 A corresponde a la Imea horizontal m = 0, que conecta la epoca 
de materia (r, m) = (—1,0) con el punto de de Sitter Pi (r, m) = (—2,0). Una posible generalizacion 
del modelo ACDM 

f{R) = {R'-Ar, (5.55) 
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Figura 5.2: Piano (r, m) para f{R) = (R^ - A)". 



con lo cual r = -c{R^ - A) ^bR^ y 

(I - c) 

m(r) = ^ '-r + b-1. (5.56) 

c 

Si la interseccion m = —1 + bc con la Imea critica es en < m ^ 1 y la pendiente esta dada por 
— 1 < (1 — c)/c < 0, entonces la epoca de materia esta conectada co Pi y el modelo es aceptable. En la 
Figura 5.2 se muestra el piano (r, m) para f{R) = {R^ — A)^. 



5.3.2. f{R) = R + aR-'' 

Este modelo esta propuesto en [66] para explicar la actual expansion acelerada. Para este modelo 
tendremos 

m(r) = _^^(i±l), (5.57) 
r 

vemos que tal modelo es independiente de a, y puesto que m{r = — 1) = el modelo satisface la condi- 
cion necesaria para la existencia del punto de materia P5. Reemplazando (5.29) para (5.57), obtenemos 
la solucion nia = o rrih = — (n + 1), que es cierto para los puntos P5 y Pg- En este caso tales puntos 
estan dados por 

Aa= (^0,-^,0, nm = l, Weff = 0, (5.58) 

3(n + l) 4n + 3 4n + 3\ Sn^ + 13n + 3 1 
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2(n + 2) 



4n + 5 



n(4n + 5) 



0, 



6n^ + 7n-l 



(5.60) 



2n + 1 ' (n + l)(2n + 1) ' (n + l)(2n + 1) 

"""■^■^ " ~ ~ 3(n + l)(2n+l) 

Esta familia de modelos se puede dividir en tres casos: (l)n < —1, (2) — 1 < n < 1 y (3)n > 1. Para el 
caso (1) tendremos que m' = n/r'^, implica que m'{r = —1) < —11 y por lo tanto la epoca de materia 
alrededor de m ~ 0^ es estable (Pi es estable tambien para 2 < n < 0). El caso con n = —2 corresponde 
al modelo de inflacion de Starobinsky [21]. 



(5.61) 



Para el caso (2) la condicion r = — 1 se satisface con P — oo y vemos que 

_ n(n + l)aP"""^ 
~ (1 -naP-"-i) ' 

se aproxima a cero cuando si a = 0. Para este caso existen oscilaciones alrededor de la epoca de materia 
y un punto final de de Sitter Pi. Para el caso (3) el punto Pq es estable y acelerado en la region (A) 
(para Weff diferente de un valor phantom) debido a la condicion m = —n — 1 < — 1. Si a > 0, m se 
aproxima a cero por el lado derecho, y entonces existen oscilaciones alrededor de la epoca de materia, 
pero el punto de aceleracion Pi es inestable. Cuando a < 0, m se aproxima a cero por el lado izquierdo 
y de nuevo el punto de materia es inestable, pues uno de los autovalores muestra una divergencia. En 
la figura 5.3 se muestra el piano (r, m) para el modelo /(P) = P + a/P, (n = 1) donde se aprecia 
claramente que NO es un modelo viable, pues no conecta el punto de la epoca de materia con el de 
expansion acelerada. En la figura 5.4 se muestra el piano (r, m) para el modelo de Starobinsky. 

f(R)=-+R 
R 



S 





Figura 5.3: Piano {r,m) para f{R) ~ R + a/R. 
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f(R) =R + aR- 



3 



2 



S 













f{R] = R-2A 




i 

.1 











-2-10 1 2 3 

r 

Figura 5.4: Piano (r, m) para f{R) = R + oiR^ , modelo de Starobinsky. 

5.3.3. f{R) = RPei^ 

En este modelo m esta dado por 

P 

m(r) = — r H — , (5.62) 
r 

En el caso exponencial puro {p = 0) se tiene m = — 1 y xs/m — >• 2:2 — )■ de modo que el punto P2 existe 
mientras que P5 no. Por otro lado vemos que la funcion m se hace cero para r = ±y^, lo que implica 
que la condicion m{r = —1) = para la existencia de un punto de materia se cumple para p = 1. Sin 
embargo, puesto que en este caso m'(r = —1) = — 2 < —1, el punto P5 es una espiral estable para m > 0. 



En el limite m — t- 0"*", el punto Pg no puede ser usado para la expansion acelerada, ademas el pun- 
to -P5 es estable. Por otra parte, cuando m(r = —2) = 3/2 para p = 1, el punto de de Sitter no es 
estable. Estos modelos entonces no tienen una secuencia de epoca de materia y expansion acelerada 
para p = 1. En la figura 5.5 se muestra el piano (r, m) para f{R) = Re^ {p = Q = !)• 



5.3.4. f (R) = RP (In aRy 

Para este modelo tenemos la relacion 

p"^ + 2pr - r{q - r + qr) 

m{r) = , (5.63) 

qr 

puesto que m(r = —1) = —{p — l)^/g, la epoca de materia existe solo para p = 1. Cuando p = 1 se 
tiene m(r = —2) = 1 — l/2q, lo que significa que Pi es estable para q > pero no lo es para q < 0. La 
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f(R) = e''R 



S 









\ 




\ fiR) = R-2A 




-\- I- 











r 

Figura 5.5: Piano {r,in) para f{R) = Re^. 

derivada de m{r) es 

r2 - 1 

m'(r) = -l + ^, (5.64) 

como m'{r = —1) = —1, el punto P5 es estable. Similarmente el punto de aceleracion Pg es estable en 
la region (C) para q > 0, pero no lo es para q < 0. Esto implica que la curva m(r) no atraviesa el 
punto Pq en la region (C), de modo que las linicas trayectorias posibles son desde el punto de materia 
P5 al punto de de Sitter Pi, (r = —2). En la figura 5.6 se muestra el piano (r, m) para f{R) = RlnR, 
{p = a = q=l). 



5.3.5. f{R) = RPei/^ 
Para este modelo se tiene que 

p + r{2 + r) 

mir) = , (5.65) 

r 

la cual es independiente de q. Tenemos que m{r = —1) = p — 1 de modo que la epoca de materia existe 
para p = I. En este caso se tiene m(r) = — (r + l)^/r > para r < 0, y puesto que m{r = —2) = 1/2 
para p = 1, el punto Pi es una espiral estable. La derivada de m(r) es m'{r) = — 1 + que implica 
entonces que m'(r = —1) = y m'(r < —1) > —1. Esto muestra que el punto P5 es un punto de silla 
mientras que Pg en la region (C) es estable. La curva m(r) satisface m(r) < — r — 1 en la region r < — 1 
para p = 1 y tambien un comportamiento asintotico m{r) — ?■ — r en el Kmite r — t- —00. Entonces, en 
principio, es posible tener una secuencia P5 — Pg (r — )• —00), pero la trayectoria desde el punto P5 
es atrapada por el punto estable de de Sitter Pi que existe en {r,m) = (—2, 1/2). En la figura 5.7 se 
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f(R) = R log(R) 



2 



S 



-1 



-3 





Pl\ 










t—r = -2 














^ 

f(R] = R-2A 








— m = -r-l 



-1 1 

r 



Figura 5.6: Piano (r, m) para /(/?) = i?lni?. En la figura se debe entender el nombre log como el logaritmo natural In. 
muestra el piano (r, m) para /(ii) = Re^l^. 

Si observamos ahora las figuras 5.2 a 5.7 encontramos que las linicas que nos pueden dar una epoca de 
materia seguida por una expansion acelerada son los modelos /(i?) = R\uRy f{R) = Re^^^. La razon 
de estos es que los demas modelos atraviesan la region delimitada por m = —r — 1 para m > — 1 hasta 
alcanzar el punto de de Sitter ubicado en la Imea r = — 2. 

Podemos entonces resumir las condiciones que debe cumplir una funcion f{R) para ser cosmologi- 
camente viable: 

1. Un modelo f{R) debe tener una epoca estandar de dominio de materia solo si satisface las con- 
diciones 

m(r) « 0^, m'{r) > — len r ~ —1, (5.66) 

2. La epoca de materia es seguida por una aceleracion tipo de Sitter {w^ff = — 1) solo si 

< m{r) < len r = —2, o m(r) = — r — 1 — t- ±00, (5.67) 

3. La epoca de materia es seguida por un atractor con w^ff > — 1 solo sim = — r — lyse cumplen 
simultaneamente 

(\/3 - l)/2 < m(r) < 1 y m'(r) < -1. (5.68) 
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Capi'tulo 5. Modelos Cosmologicos en Teon'as f{R) 




En la siguiente seccion estudiaremos el problema de las distancias cosmologicas en teon'as f{R) estu- 
diando para esto la Ecuacion de Desvio Geodesico en estas teorias. 
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Capitulo 



6 



Distancias Cosmologicas en Teorias f{R) 



onsideramos ahora el problema de las distancias cosmologicas en teorias f{R)- Este estudio es de 



gran importancia, pues permite restringir los modelos f{R) mediante observaciones, en particular 
por ejemplo con la medicion de distancias en supernovas de tipo lA [3]. Si bien, las ecuaciones obtenidas 
en el Capitulo 5 nos permiten en principio determinar las distancias, usando tambien las expresiones 
dadas en el Capitulo 3, debemos notar sin embargo que tal proceso conlleva una complejidad inherente 
en cuanto la misma forma de la funcion f{R)- 

La expresion (5.4) nos da entonces una relacion para el parametro de Hubble H en funcion de la 
densidad de energfa p, la presion p y fuertemente en el escalar de Ricci R, por otro lado, sin embargo, 
de la metrica de Robertson- Walker se tiene una relacion para R en funcion de H y H, y ademas una 
relacion diferencial para R en funcion de la traza del tensor energfa-momentum, a traves de la traza de 
las ecuaciones de campo (4.36). De esta forma es necesario o bien fijar la forma para H{z), puede ser por 
observaciones o suponiendo un ansatz, o bien fijar la forma de la funcion f{R) para poder determinar 
H{z) partiendo de las ecuaciones modificadas. 

En el formalismo de Palatini, por ejemplo, las distancias cosmologicas pueden ser determinadas y 
comparadas con datos observacionales, fijando la forma de la funcion f{R) [81]. Cabe decir que la for- 
ma de las ecuaciones de campo en el formalismo de Palatini permiten en principio determinar la forma 
para H{z) de una manera mas simple que en el caso del formalismo metrico, debido a que en estas 
ecuaciones los operadores □ y V^V^ no aparecen de forma explicita en las ecuaciones. 

Un metodo interesante para atacar el problema de las distancias cosmologicas en teorias f{R) es pre- 
sentado en [82] a traves de la denominada cosmografia. El punto clave en la cosmograffa es la definicion 
de nuevas parametros cosmologicos que contienen derivadas del factor de escala hasta de quinto orden. 
Ya conocemos uno de ellos, el parametro de Hubble dado por H{t) = ^, ademas definimos [83] 




a 



Parametro de desaceleracion 
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Parametro de jerk (6.1b) 
Parametro de snap (6.1c) 

Parametro de lerk (6. Id) 

con lo cual el factor de escala a{t) se puede escribir como una serie en terminos de t. Estos parametros 
nos permiten ademas obtener las expresiones para las distancias (diametral angular y de luminosidad) 
como series en funcion del redshift [83]. Aunque esta aproximacion tiene algunas ventajas (es cosmologia 
sin ecuaciones de campo) conlleva al problema de que nos da tan solo una aproximacion a los valores 
para las distancias, dependiendo los ordenes usados en la expansion. El metodo seguido entonces en 
[82], extiende tal metodo para la determinacion de las distancias en teorias f{R) fijando la forma de 
f{R) como una expansion en serie de Taylor. 



jit) 

S{t): 

m 



+ 



a 



+ 



+ 



Nuestro objetivo en este capitulo es utilizar al EDG extendida a teon'as f{R) para la determinacion de 
distancias de una forma alternativa. 

6.1. Ecuacion de Desvio Geodesico en Teorias f{R) 

Como vimos en las secciones 2.3 y 3.3 el estudio de la Ecuacion de Desvfo Geodesico (EDG) nos per- 
mite caracterizar las propiedades mas importantes de un espacio-tiempo, y lo mas importante tratar 
el problema de las distancias cosmologicas de una manera muy general. Nuestro objetivo es entonces 
extender los resultados de la EDG para las teorias f{R) en el formalismo metrico siguiendo para esto 
[84]; en el formalismo de Palatini la EDG se ha estudiado en [85]. 

Nuestro punto de partida es entonces son las expresiones para el tensor de Ricci y el escalar de Ricci 
escritos a partir de las ecuaciones de campo y su traza respectivamente (4.35) y (4.36) 



Rn 



1 



f'iR) 



R 



f'{R) 



KT + 2f{R)-3nf{R) 



(6.2) 
(6.3) 



Usando estas expresiones podemos escribir el tensor de Riemann (3.34) como 



2f>{R) 



7ai 



6f{R) 



Tsagp-y) + f{R) {ga-ygsp - gasgj/s) 



kT + 2f{R) + 3n/'(i?) ) {ga-ygsp - gasg^p) , 



(6.4) 
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donde definimos el operador 



(6.5) 



Entonces subiendo el primer mdice en el tensor de Riemann 
1 



-«/37<5 - '-^/37<5 + 2f'{R) 



1 



6/'(i?) 



kT + 2f{R) + 3n/'(i?) {6^gs/, - 6fg^is), 



(6.6) 



y contrayendo con 

y/3^7y<5 el lado derecho de la EDG (2.78) nos queda 



2f'{R) 



+ {6^Vsp - + gpsV^^ - gp^V,^)nR) 



- (y^T + 2f{R) + 3n/'(i?)) {6'^gsp - 5'ig,p)V^r,'^V' , (6.7) 

Esta ecuacion nos permitira entonces encontrar la EDG para cualquier espacio-tiempo y para cualquier 
forma del tensor de energia-momentum en teorias f{R)- 

En la siguiente seccion mostraremos en detalle los pasos para encontrar la EDG en teorias f{R) usando 
la metrica de Robertson- Walker, nuestro proposito tambien es el de comparar los resultados con los de 
RG para el caso Ifmite f{R) = R — 2 A. 

6.1.1. Ecuacion de Desvio Geodesico para el Universe FLRW 

Usando entonces la metrica de Robertson- Walker (3.1) y el tensor de energia-momentum para un fluido 
perfecto (3.2) tendremos 



Rn 



1 



fiR) 

R = 



fiR) 

k{p + p)UaUi3 + ( H — ) gal3 + Vapf{R) 



1 



fiR) 



K{3p-p) + 2f{R)-3nf'iR) 



(6.8) 
(6.9) 



con estas expresiones podemos escribir el tensor de Riemann como 
1 



R. 



n{p + p) {usUi3gaf - U^UpgaS + U^UagpS - USUagpf) 



2/'(i?) 

+ [ I^P+'y ~^ + ^f'iR) ) {ga-ygSfS - gaSg-yp) + {ga-yT^Sp " + 5/35^7^ " g^-yl^Sa) f (R) 



3.10) 
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donde de nuevo hemos usado que para este espacio-tiempo el tensor de Weyl Ca/S'yS es nulo. Ahora, si 
el campo vectorial V" esta normalizado, tendremos que V^Va = e, y por lo tanto 



1 



subiendo el primer mdice en el tensor de Riemann y multiplicando por rf' 
1 



3.11) 



R%sV^r^^V' 



2f'{R) 



+ (^Kp+^ + ^-^+Uf\R)yer-V''{V-,r,^)) + [{5^V,p-5'^V,p+g 



con E = —VaU'^-.Tjau"' = rjaV"' = [35] esta expresion se reduce a 
1 



5^7 



3.12) 



R'^p^sVri^V' 



2f'{R) 



+ 



1 



r]\ (6.13) 



2f{R) 

Usando de nuevo la expresion (6.14) 

r = q\^ + (^Y + \\ =6\h + 2H' + \\, (6.14) 
a \a J a^J [ 

Usando ahora los resultados del Apendice C tenemos que los linicos operadores que son diferentes 
de cero son 



nf{R) = -dlf'{R) - SHdof'iR), 

= -f"{R)R - f"'{R)R? - mf"{R)R, 

Poo = -^HdofiR), 
= -mf"{R)R, 

= 2Hgijdof'{R) + g:jdlf{R), 

= 2Hgijf"{R)R + gij{r\R)R + f"{R)R?), 



3.15) 



3.16) 



3.17) 



con gij las componentes espaciales de la metrica de Robertson- Walker y R = d^R. Con estos resultados 
podemos obtener la contribucion total de los operadores en (6.13), ver Apendice D 

{6'^V,p - dfV^p + g^sV^^ - gp^V,^)f{R)V^V''n^ = e{f>Hf"{R)R + f"{R)R + /"(i?)^^)^, (6.18) 
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entonces la expresion para R'^^^V^rj'^V^ se reduce a 



2f'{R) 



, fiR) 



2Hf"{R)R 



V 



3.19) 



la cual es la generalizacion de la ecuacion de Pirani. En el caso particular f{R) = i? — 2A la anterior 
expresion se reduce a (3.43). 



Entonces podemos escribir la EDG en gravedad f{R) a partir de la expresion (2.78) 



1 

'2f'{R) 



K{p + p)E^ + e( + ^ + + 2HnR)R 



3.20) 



como esperabamos la EDG induce solo un cambio en la magnitud del vector desviacion ij", lo cual 
ocurre tambien en RG. Este resultado era de esperarse debido a la forma de la metrica, que describe 
un universo homogeneo e isotropico. Para universos anisotropicos, como por ejemplo los de Bianchi I 
la EDG induce tambien un cambio en la direccion del vector desviacion, como se muestra en [39], [86]. 



6.1.2. EDG para Observadores Fundamentales 

En este caso tenemos V"' como la cuadri-velocidad del fluido n". El parametro afm v coincide con el 
tiempo propio de los observadores fundamentales u = t. Puesto que tenemos geodesicas temporales 
entonces e = — 1 y tambien E = 1, entonces a partir de (6.19) 



1 



2f'{R) 



2kp f{R) 



2Hf"{R)R 



7] 



si el vector desviacion es ija = ica, la isotropi'a implica 



0, 



Dt 

y 

~dW ~ ~d^^ ' 

usando este resultado en la EDG (2.78) con (6.21) tenemos 

df^ 



1 



2f'{R) 



2kp f{R) 



2Hf"{R)R 



en el caso i = a(t) tendremos 



a _ 1 
a " i%R) 



3.21) 



3.22) 



3.23) 



3.24) 



3.25) 



Esta ecuacion puede ser obtenida como un caso particular de la ecuacion de Raychaudhuri generalizada 
dada en [87]. Es posible obtener como ejemplo, la forma estandar de las ecuaciones modificadas de 
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campo a partir de la ecuacion de Raychaudhuri. Primero restando la expresion R/6 en (6.25) con R 
dado en (6.14) tendremos 



a a 
a a 



1 



k 



a2 f{R) 



^+HriR)R-^ 



k_ 
k 



1 



Kp f{R) 



f'iR) 
1 

1 



Hf"{R)R + 



3 6 

- {Rf'{R)-f{R)) 



R 

Rf'jR) 

6 



3Hf"{R)R 



3/'(i?) 



3.26) 



que corresponde a la ecuacion (5.4). Para obtener la segunda ecuacion sumamos a ambos lados de la 
ecuacion (6.25) la expresion R/6, y usamos tambien la forma para la traza de las ecuaciones de campo 
(6.3) con el operador □/'(i?) dado en (D.2) 



a 


d 




2 


k 




+ - + 


(fly 


) - 




a 


a 






^ ~ 




d 




2 


k 




2- + 




) - 






a 


V a , 







3 



+ 



R 

6' 



6 3 3 ^ 3 3 V ; > 



f"'{R)R^ -2Hf"{R)R + 



RF{R) 
6 



2H + 3H'^+ ^ 



a2 /'(i?) 



2F + 3H^ + 



A; 



-Kp - 2Hf"{R)R - (/(-^) Rf'iR)) _ f"(^R)ji _ /'"(i?)i?2 
Kp + 2Hf"{R)R + (/(-^) - Rf'iR)) ^ ^//(^)^ ^ /"'(i?)i?' 



3.27) 



que corresponde a la segunda ecuacion de Friedmann modificada (5.5). 
6.1.3. EDG para Campos Vectoriales Nulos 

Consideramos ahora la EDG para campos vectoriales nulos dirigidos al pasado. En este caso tenemos 
ya _ j^a^ kak" = 0, entonccs la ecuacion (6.19) se reduce a 



1 



K{p + p)E\", 



2f'{R) 



3.28) 



expresion que puede ser interpretada como el Enfocamiento de Ricci en teorias f{R)- Escribiendo tj" = 
Co-e" = 1, eau" = Cak" = y escogiendo una base paralelamente propagada ^jy^ = k^Vj3e°' = 0, 
la EDG (6.20) nos queda 

(frj 



1 



-k{p + p)E'^ t]. 



3.29) 



V'{R) 

En el caso de RG discutido en [35] , todas las familias de geodesicas nulas dirigidas al pasado experimen- 
tan enfocamiento, siempre y cuando se cumpla n{p + p) > 0, y para un fluido con ecuacion de estado 
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p = —p (constante cosmologica) no hay influencia en el enfocamiento. A partir de (6.29) la condicion 
para enfocamiento en teorias f{R) es 



f'iR) 



> 0. 



3.30) 



Una condicion similar sobre la funcion f{R) fue establecida para evitar entre otras la aparicion de 
fantasmas ^ [11], [21], en la cual f'{R) > 0, ver tambien Seccion 5.1. Queremos ahora escribir la ecuacion 
(6.29) en funcion del parametro de redshift z. Siguiendo la Seccion 3.3.2 podemos escribir entonces 



dv _ 1 

Tz ~ EqH{1 + zY' 

1 



EqH{1 + z) 



1 ,dH 



J TT 

expresamos nuevamente como 



dH dv dt dH 



1 



dH 



dz dz dv dt H{1 + z) dt 



De la definicion para H tendremos 



H 



dH da a 



dt dt a a 
and usando la ecuacion de Raychaudhuri (6.25) 



H' 



H 



1 



m 

6 



Hf"{R)R 



Kp 



H' 



entonces 



d'^v 



dz^ EqH{1 + zf 



el operador 4-7 toma la forma 



dv 



dv^ 



d^rj 3 

+ 



(1 + z) 



dr] 
dz 



y la EDG (6.29) se reduce a 



d^7] 

d^ ' (l + z) 



+ 



dz 2m{l + zYf'{R)'^ ' 



3.31) 
3.32) 

3.33) 

3.34) 

3.35) 
3.36) 



(6.37) 



(6.38) 



Ahora, la densidad de energi'a p y la presion p considerando las contribuciones de materia y radiacion, 
se pueden escribir como 



KP = 2,Hlnm0{l + zf + 2,Hlnro{l + z)^ Kp = i?o^Or-o(l + zf, 



3.39) 



^En donde fantasmas se refieren a estados con norma negativa o campos con signo negativo en el termino cinetico. 
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donde hemos usado pm = y Pr = ^Pr- Podemos ahora escribir la EDG en una forma mas compacta 



p^+r{H,R,z)^ + QiH, R, z)^ = 0, 
az'^ az 



con 



V{H,R,z) 



(1 + z){nmo{l + z)^ + 0,o(l + z)^ + f'{R)nko{l + zf + VLde) 



QiH,R,z) 



3nm0{l + Z) + 4nro{l + Z) 



2{nmo{l + z)^ + 0,o(l + z)^ + f'{R)nko{l + zf + ^de) ■ 



y H dado por las ecuaciones de campo modificadas (5.4) 



1 



f'{R) 



H"^ = Hi 



Hln^oil + zf + HlQMl + zf + f^^^^ - Hf"{R)R 



{nmoil + zf + r?,o(l + zf + Qde) + ^k{l + zf 



f'{R) 



donde 



DE 



Hi 



6 



Hf"{R)R 



que es equivalente a la expresion (5.23) y 



3.40) 



3.41) 



3.42) 



3.43) 
3.44) 



n 



kO 



H'2„2 • 



3.45) 



Para poder resolver la ecuacion (6.40) es necesario escribir R y H en funcion del redshift. definimos 
entonces 

d _ dz da d _ ^ d 
dt da dt dz dz ' 

y entonces para el escalar de Ricci tendremos [7] 

2 ; 



R = 6 
= 6 
= 6 



a a 
- + (- 
a \a 



+ 



2H^ + H + 



2H^-{l + z)H^ + k{l + z) 



Ahora, para tener H = H{z) o bien sea fijar de entrada la forma H{z) o fijar una forma especifica de 
la funcion f{R). Este punto ha sido discutido en detalle en [10] y el metodo de fijar H{z) para hallar 
la forma de la funcion f{R) por observaciones en 



En el caso particular f{R) = R — 2A tenemos f'{R) = 1, f"{R) = 0. la expresion para VtoE se 
reduce a 

'{R-R + 2kf 



DE 



1 



6 



A 

3Hl 



3.47) 
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y por lo tanto, como se habia anticipado en el Capitulo 5 el parametro ^de generaliza el parametro de 
energia oscura La ecuacion de Friedmann modificada (5.4) se convierte en la expresion conocida en 
RG (3.23) 

= Hi[nmoii + zf + nrQii + z)^ + nA + nk{i + zf], (6.48) 

las expresiones 7^, y Q se reducen a 

. ^ 417,o(l + z)4 + (7/2)1^^0(1 + zf + 317fco(l + zf + 20a 
(1 + z)(l7,o(l + ^)^ + f^mo(l + 2)3 + f^fco(l + ^)2 + ^^a) ' 

^( ^ ^ 21^,o(l + z)^ + (3/2)17^0(1 + z) 

y la EDG para campos vectoriales nulos es 

A 417^o(l + z)^ + (7/2)17^o(l + zf + SQkoil + ^)^ + 217a dr] 

dz^ ^ (1 + Z)(0^o(l + Z^ + 17™o(l + Z)3 + Ofco(l + 2)2 + Qa) dz 

217,0(1 + 2)2 + (3/2)17^0(1 + 2) ^ 

^17,o(l+^)4 + 17™o(l + 2)3 + 17fco(l + 2)2 + 17A'' • ^ - ^ 

que coincide con la expresion (3.76). Nuevamente la relacion de Mattig es obtenida en el caso 17a = 
y escribiendo 17^0 = 1 — 17mO — ^ro 

A 6 + 17^o(l + 7z) + a.o(l + 8z + 4z^) dri 317,^o + 417,o(l + 2) ^ ^ 

dz^ 2(1 + z)(l + 17^02 + a,oz(2 + z)) dz 2(1 + z)(l + 17^02 + 17,02(2 + z)) ' ^ ' ' 

de este modo la ecuacion (6.40) nos da una generalizacion a la relacion de Mattig en teorias f{R)- 

Para relacionar ahora la magnitud del vector desviacion rj con la distancia diametral angular, con- 
sideramos nuevamente que la magnitud del vector desviacion rj se relaciona con el area propia dA de 
una fuente a redshift z por dri oc VdA, y de esta expresion, la definicion de la distancia diametral 
angular dA puede ser escrita como 

con dl7 el angulo solido. De modo que la forma de la relacion de Mattig generalizada es la misma 
cambiando r] por d^ (el factor de proporcionalidad se puede descartar de la ecuacion diferencial) , de 
modo que 

^2 ^m) 417„o(l + zf + 417,o(l + 2)4 + 3/'(i?)17fco(l + zf + AQde - ^ff- dd{^^^ 



dz^ {i + z){nmo{i + z)^ + nro{i + zY + f'{R)nko{i + zy + nDE) dz 

317^o(l + 2) + 417,o(l + z)2 ^y(R) 

2(17^0(1 + 2)3 + 17,o(l + 2)4 + /'(i?)17fco(l + 2)2 + ^de) ^ 



^Cabe aclarar que la generalizacion es para contribuciones netamente geometricas. 
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f( 

donde denotamos la distancia diametral angular por , para enfatizar que alguna solucion de la 
ecuacion anterior necesita una forma especifica de la funcion f{R), o bien la funcion H{z). Esta ecuacion 
satisface las condiciones iniciales (para z > zq) 

d^}^\z,ZQ 



dd 



f{R) 
A 



dz 



-{z,Zq] 



= 0, (6.55) 

Z=ZQ 

H{zq)(1 + zq) 



con H{zq) siendo la ecuacion de Friedmann modificada (5.4) evaluada en z = zq. 

6.2. Relaciones entre Distancias Cosmologicas en Teorias f{R) 

En el Capitulo 3 mostramos las diferentes relaciones entre las distancias cosmologicas en RG. Nos 
preguntamos entonces si tales relaciones se mantienen validas en las teorias f{R)- Si bien tales relaciones 
entre las distancias las mostramos definiendo en si cada una de ellas, existe un teorema denominado 
teorema de Reciprocidad mostrado por Etherington en (1933) [89] y que nos dice 

Teorema 6.2.1 (Teorema de Reciprocidad (Etherington)). La distancia area observador ro y la dis- 
tancia area galaxia re son iguales salvo un factor de redshift cosmologico 

rl = rl{l + zf. (6.57) 

El punto fundamental de este teorema es el hecho de que varias propiedades geometricas son invariantes 
cuando los roles de fuente y observador se intercambian. La distancia re se relaciona con la distancia 
de Luminosidad por la relacion 

dL = rG{l + z), (6.58) 

y por lo tanto 

di = ro(l + z)2, (6.59) 

y entonces usando que dA = r^dO. podemos relacionar la distancia area observador con la distancia 
diametral angular, para obtener 

dL = dA{l + zf. (6.60) 

Una prueba de este teorema de reciprocidad esta dado en [35], [86] a partir de la primera integral en la 
EDG. Mostraremos a continuacion muy brevemente estos resultados. 

6.2.1. Primera Integral en la EDG 

La primera integral nos relaciona basicamente los vectores desviacion ry", r/2 de dos geodesicas que 
comparten el mismo campo vectorial geodesico V^. Es decir, ambas cumplen a partir de (2.78) 



R%5VW,V\ (6.61) 



^ = (6.62) 
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Multiplicando (6.61) por r]2a Y (6.62) por r]ia J teniendo en cuenta que podemos subir y bajar los 
indices que se suman 



= -V2aR''p,sV^r^JV' = -v^RafS^sV^vjV, (6.63) 
^1°^ = -ViaR'-p.sV^vlV' = -7ifR^^,sV^r,JV', (6.64) 



teniendo en cuenta que las derivadas covariantes cumplen la regla de Leibnitz podemos escribir 

Dv'^ Dv \ Dv / Dv Dv 



Via-^ = 7^ - ^„ ^„ . (6.66) 



Dv"^ Dv \ ' Dv J Dv Dv 

restando estas expresiones 



D'V^ irr]l_D_f Dr^_ Dr]f\ f Drj2a Di^f Drj.p, Drj^ 



D 


( Drfi 


Dt]? 


Wv 


l^^" Dv 


Dv 


D 


( Dt]^ 


Dn? 


Dv 


['^"^ Dv 


Dv 



(6.67) 

Por otro lado, restando (6.63) de (6.64) y teniendo en cuenta las propiedades de simetria del tensor de 
Riemann tendremos 

= Ra^-ysVf'V'ivWl) - R^SapV^ViVlVl) 

= Ra^-ysVf'v'ivWi) - Ra^SjsV'V^irjJrj^) 

= 0, (6.68) 



de modo que 

Dv V Dv Dv 
con lo cual finalmente obtenemos la primera integral de la EDG 



3.70) 



Dr]^ Dt]? 

Vla-f:: ma-f^ = COnst. 

Dv Dv 



Esta ecuacion es valida para cualquier modelo cosmologico que supone valida la Relatividad General. 
Esta primera integral se puede investigar para el caso de rayos nulos. Consideremos geodesicas nulas 
que diver gen de una fuente S y llegan a un observador O, con vector desviacion rji y rayos nulos que 
divergen de la fuente S y llegan al observador O, con vector desviacion r]2. La primera integral esta 
expresada en terminos del parametro affn v, pero expresemosla ahora en terminos de la magnitud i 
definida por di := Sdr. Para encontrar la relacion entre dr y dv, consideremos una geodesica radial de 
la siguiente manera 

y":=Sn" + Pe", (6.71) 
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donde e° = S" ^(9r)° = -S* ^(^" , eo-e" = 1, eau" = 0. Para una geodesica radial nula tenemos ademas 

^ = V' = E = -{V^u^ = Eo^, (6.72) 
dr 1 P{S) E{S) 



De esta forma 



podemos escribir entonces 



M = Sdr = Eo[^] du, (6.74) 



dr]i didr]i So drji ^ . .drji 

-^ = m + ^)^. (6.76) 
Teniendo en cuenta que para el observador z = tenemos 

??2|o-^|o = r/i|s-^|5(l + ^;), (6.77) 

donde los terminos ^ son los angulos subtendidos por los pares de rayos nulos y corresponden a los 
vectores desviacion. Expresado en terminos de las distancias area y luminosidad, definidas por 

Vi\s ■■= ro-^\o, (6.78) 

V2\o ■= rc-j^ls- (6.79) 

encontramos el famoso teorema de reciprocidad para los modelos de FLRW 

rG = ro{l + z). (6.80) 

Ahora, puesto que usamos como metrica la de Robertson- Walker y las expresion para el redshift siguen 
siendo las mismas que en RG, las relaciones entre las distancias cosmologicas se preservan 
tambien en las teorias f{R)- Asi, tendremos entonces para la distancia de Luminosidad di en 
teon'as f{R) que 

(6.81) 



El proceso entonces para determinar las distancias cosmologicas en teorias f{R) se puede describir como 

1. Fijar una forma para H(z) o bien una forma especifica de la funcion f{R)- 

2. Encontrar a partir de las ecuaciones de Friedmann modificadas las expresiones para H{z) y R{z), 
en general a partir de un calculo numerico 

3. Determinar Qde (6.44) en funcion de z. 

4. Resolver la ecuacion (6.54) para d } con las condiciones iniciales (6.55) y (6.56). 
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5. Obtener a partir de esta expresion las demas distancias cosmologicas. 

Si bien, nuestro metodo para la determinacion de distancias a partir de la EDG para campos vectoriales 
nulos en teorias f{R) necesita una forma especifica para la funcion, es mia herramienta elegante y que 
nos brinda la posibilidad de comparar directamente con resultados observacionales, siendo asi una buena 
herramienta para restringir modelos f{R)- 



6.3. Ecuacion de tipo Dyer-Roeder en Teorias f{R) 

Finalmente daremos una importante relacion que nos permitira estudiar las distancias cosmologicas 
tambien en universos con inhomogeneidades. La ecuacion de Dyer-Roeder nos da una ecuacion diferen- 
cial para la distancia diametral angular dA en funcion del redshift z introduciendo un parametro que 
nos permite introducir inhomogeneidades en la densidad de materia [90], [91]. La forma estandar de la 
ecuacion de Dyer-Roeder en RG esta dada por [40], [92] 

(1 + ^)'-^(^)^ + (1 + ^)^(^)^ + '^i^)dA = 0, (6.82) 

con 

T{z) = H'^{z), (6.83) 
dH 

g{z) = {l + z)H{z)— + 2H\z), (6.84) 

nz) = ^^n^o{l + zf, (6.85) 

siendo a{z) el parametro de suavidad, que nos da el caracter de las inhomogeneidades en la densidad 
de energia. 

Con el fin de obtener una ecuacion de tipo Dyer-Roeder en teorias f{R) seguimos [93]. Primero de- 
bemos notar que los terminos que contienen las derivadas de en la ecuacion (6.54) vienen de la 

transformacion ^ — > ^ y el termino que acompaiia a d^^^^ viene del enfocamiento de Ricci (6.28). 
Entonces siguiendo [40], [91], introducimos la fraccion de masa q (parametro de suavidad), y entonces 
reemplazamos solamente en el enfocamiento de Ricci p — > ap. Siguiendo estos argumentos en (6.40), 
y considerando el caso QrO = tendremos 

+ d? + + (o^o(l + ^)=^ + /'(i^)^^fco(l + ^)2 + ^^DE) d~z 

, 3a{z)^rn0{l + ^f{R), DR ^ q gg. 

2{am0{l+z)^ + f'{R)Qk0{l + z)^ + nDE) ^ • V - ^ 

donde denotamos la distancia de Dyer-Roeder en gravedad f{R) por d^J^^^' y ademas a diferencia 
de la expresion (6.54), cada uno de los terminos esta dividido sobre {1 + z)"^. Esta ecuacion generalizada 
de Dyer-Roeder tambien satisface las condiciones (6.55) y (6.56), y como es de esperar, se reduce a la 
forma estandar en RG para el caso particular f{R) = R — 2A. 
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Capitulo 



7 



Conclusiones y Perspectivas 



omo vimos a lo largo de este trabajo, las teorias de gravedad modificada f{R) proveen una intere- 



sante alternativa para atacar el problema de la actual expansion acelerada del universo. Su interes 
creciente en los ultimos 15 afios radica en parte por ser una de las generalizaciones mas simples de la 
Relatividad General, consistiendo de una funcion arbitraria del escalar de Ricci, en el lagrangiano que 
describe la gravitacion. 

En el Capitulo 4 obtuvimos las ecuaciones de campo en el formalismo metrico de las teorias f{R) 
destacando la importancia de los terminos de frontera en la accion modificada. Las ecuaciones de cam- 
po obtenidas nos muestran claramente que son de cuarto orden la metrica, y tambien como pueden 
interpretarse los terminos adicionales, como la contribucion a un tensor de energia- momentum efectivo 
compuesto de componentes netamente geometricas. Discutimos tambien las ecuaciones de campo en los 
formalismos de Palatini y metrico-afm, y la equivalencia de las teorias f{R) con las de Brans-Dicke. 

Los modelos cosmologicos fueron estudiados en el Capitulo 5, partiendo para esto de las ecuaciones 
de campo modificadas y de la validez del principio cosmologico, lo que nos permitio usar como espacio- 
tiempo el de Robertson- Walker. El estudio de las ecuaciones de Friedmann modificadas modeladas como 
un sistema dinamico autonomo, permite encontrar criterios para obtener modelos f{R) cosmologica- 
mente viables, y obteniendo ademas criterios para las funciones que conecten una epoca de dominio de 
materia, con una de expansion acelerada. 

El estudio de la Ecuacion de Desvfo Geodesico EDG en las teorias f{R) nos permitio tratar el pro- 
blema de la medicion de distancias, realizado en el Capitulo 6. Obtuvimos a partir de la EDG la 
generalizacion para la ecuacion de Pirani, y en el caso de observadores fundamentales las ecuaciones 
de Friedmann modificadas escritas en la forma estandar. El caso de la EDG para campos vectoriales 
nulos permitio encontrar expresiones generales para el enfocamiento de Ricci, la relacion de Mattig, 
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y la ecuacion general que determina la distancia diametral angular (Ia en teorias f{R), lo que permi- 
te determinar las distancias cosmologicas (en particular la distancia de luminosidad y la distancia 
comovil x{z))- Finalmente obtuvimos una expresion para la ecuacion de tipo Dyer-Roeder en teorias 
f{R), abriendo la posibilidad del estudio de distancias cosmologicas en universos con inhomogeneidades. 

Si bien las teorias f{R) sufren algunos problemas en cuanto a su viabilidad en pruebas del sistema 
solar [114], tiene pruebas a grandes escalas que las hacen ver como una posible explicacion a la actual 
expansion acelerada. Ultimamente se ban hecho extensiones de trabajos en RG a teorias f{R) (des- 
de lentes gravitacionales [115], agujeros negros [116], estabilidad para formacion estelar [117], Ifmites 
Newtonianos [118], etc.), lo que permite encontrar varias formas para la prueba observacional de las 
teorias. 
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Apendice 



A 



Evaluacion del termino 



Hemos calculado ya la variacion dV^^ 

"^r^a = ^Sg'''^ [di^g^a + daQ^p - d^gpa] + ^5'^'^ [dpiSg^a) + da{5g^ii) - d^{6gpa)\ ^ (A-1) 

escribiendo las derivadas parciales de las variaciones en la metrica con la expresion para la derivada 
covariante 

^-t^Qap = d^Sgap - V^a^g^i^ - V^pdgau, (A. 2) 

y usando tambien que estamos en una variedad libre de torsion, i.e., que el simbolo de Christoffel es 
simetrico F^^ = F"^, podemos escribir 

+ ^9'''^[^l3{Sgfa) + Va(<5ff7/3) - V^(%„) + T^Jg^x + ^afi^gX'y] , 

= \dfi9^ot + dag^fi - d^g^c\ + g^^V^pJg^x + ^g'^'^ \y fii^a-to) + ^ ai^a-i^) - Vt,(%q)] , 

(A.3) 

usando la ecuacion (2.15) en el segundo termino 

+ ^g""^ i^fsiSg-ya) + ^a{5g^l3) - V^{5gi3a)\ , 

= Sg'^-gxuT)^^ - <55'^'^<5^5A.r^„ + ^5^^ pi^g^o) + Va{Sg.,p) - V^(%«)] , 

= 55""5A.r^a - Sg'^'^gXuT^pa + \g''^ [V/3(5<77a) + V„(5<7^;3) - V^(%a)] • (A.4) 
Entonces tenemos 

<5r^a = ^5"^ [V/3(55a7) + V„(%^) - V^(%„)] , (A.5) 
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y similarmente 

5Tl^ = \g^^[V^{5g,,)\. (A.6) 

Sin embargo, es conveniente expresar el resultado previo en funcion de las variaciones dg'^^ , de nuevo 
usando (2.15) 

-^r^a = \9''^['^p{-9c.^.g^u59n+'^o.{-9p^.g^u5gn " ^{-9p>.9c.u5gn], 
= -^5"^ [9a^9^uVp{S9n + g^t.9^uVa{6gn - 9p^.9auV^{Sgn] , 
= -\ K9a^Np{59n + K9p^ya[S9n - 9pt.9au9''^V^{59n] , 
= -\ [Sa^V/^Ct^^"^) + 9M'^a{6g''^) - gp^.g.uV'iSgn] , (A.7) 
donde usamos = g^'^'V^. De forma analoga 

5Tl^ = -^g^,V^i5gn- (A.8) 
Calculamos ahora el termino g^^iSV^^) — g°"^{6ra-y) 

5'^'5(5r^j _ g-'^{Srl^) = - l[[g^''g^,Vp{6g'^^)+g''^g^,V^{6g'^^) - g^^^g^^^gauV^Sgn] 

-[g''''g^.uVa{Sgn]), 
= - ^(^[S^Vf^iSg'^'') + S^V^iSg'^'') - S'^^g^.V^idgn] 

- [gf^ug^'^VaiSgn]), 

= - ^ (2V^(<55'^T) - 2g,,V^{6g^n) , (A.9) 

tenemos entonces, 

5"^(<^r^j - ^"'^(.^r^:^) = ff^.v^i^^,?^-) - v^(<55'^^). (A.io) 
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Apendice 



B 



Integrales con Mr y N 



G 



(B.l) 



dV 



Tomando la derivada covariante en M^- 

=VAnR))g^^Vr{5g'''') + nR)g^(sn{6g^^) - VA5g''^)gc.pV rifiR)) 
-6g^Pg^pU{f{R)), 

=nR)go.pU{5g"P) - Sg^^ g^pUifiR)). (B.2) 
en donde hemos usado la compatibilidad metrica gap = 0, integrando esta expresion 

/ d^x^gWMr= I d^x^gnR)gapU{5g^^)- [ d^x ^5g^^ gapOif (R)), (B.3) 
Jv Jv Jv 

usando el teorema de Gauss-Stokes (2.26), la primera integral puede escribirse como un termino de 
frontera 

d%e^\n^Mr= [ d^x^f'iR)gapn{6g''^))- [ d^x ^Sg'^^ga^nifiR)), (B.4) 
Jv Jv 

podemos escribir entonces 

/ d^x^f'iR)gai3a{6g''^)= [ d^x ^gdg^^^ g^pU{f {R)) + i d^yev^n^M,. (B.5) 
Jv Jv JdV 

De forma similar, tomando la derivada covariante en A^"^ 
V^iV- = V4/(i?)V-,(55"^)) - V^Sg'^^V^nR))), 

= VM'{R))^i{5f^) + nR)y.yA5g^^) - VA5g^^)V.,{f{R)) - Sg^^V^V^fiR)), 
= f'{R)V^V^{6g'^(') - Sg'^^V^V^WiR)), (B.6) 

integrando 

I d^x^V^N-= [ d^x^gf'{R)V„Vp{5g''P)- [ d^x ^Sg^^V^VpifiR)), (B.7) 
Jv Jv Jv 

usando de nuevo el teorema de Gauss-Stokes podemos escribir 

/ d^xV=9/'(^)V.V;3(^5"^)= / d^x^g5g''^V„Vp{nR))+ i d^ye^\n^N\ (B.8) 
Jv Jv JdV 
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Apendice 



C 



Ecuaciones de Priedmann Modificadas 



Primero escribimos los simbolos de Christoffel no nulos para la metrica (3.1) 
donde definimos la parte espacial de la metrica por 



Qij = diagi Y 



1 



„2 2 „■ 2 , 



a 



r , r sm 



(C.l) 



(C.2) 



(C.3) 



Aunque hay mas simbolos de Christoffel no nulos, una propiedad de los espacios maximalmente sime- 
tricos nos evitara ese calculo. Usemos ahora la definicion para las componentes del tensor de Ricci 



Rap — So-r^Q, — dpV^^ + r^^r^^ — r^^r^^,, 



con lo cual 



00 



3-, Rqj = 0, Rij = [ad + 2a^ + 2k)gij, 



La curvatura total esta dada por 

R = g^^Roo + 4 ^^^R = -3° (-1) + 4 (3aa + Gd^ + = 6 
Primero vamos a evaluar los terminos ^ fif {R) y n\f'{R) 

{VaVp)nR)=Va{dpf'{R)), 

= dadpf[R)-Vlpd^f'{R), 

y 

UnR)=g''^{V,V^)nR), 

= g'''{dpd,f{R)-Tl,d,f'{R)), 
= gP''dpdJ'{R)-gP'^Tl,d^nR). 



a I a\ k 



a \ a 



(C.4) 
(C.5) 

(C.6) 



(C.7) 



(C.8) 
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Entonces la componente temporal-temporal de las ecuaciones de campo (4.35) es 

nR)Roo - \f{R)gm - [VoVo - gmW{R) = i^Too, 

- S^iR) + ^ - [VoVo + a]f'{R) = Kp, 
el termino VqVo es simplemente 



(C.9) 
(C.IO) 

(C.ll) 



(VoVo)/'(i?) = dif'iR), 

y para el termino □/'(/?) tendremos 

afiR) = g'^'dif'iR) + g'^d.d.f'iR) - g^^Tl^dofiR) - g^^Tl^d^fiR) - g^^rl,dof'{R), 



-d'^f{R)+ii3^^R) 



1 _ 1 

'ah)-^ ^dof{R) 



1 



1 — kr^ 



iaay^dofiR) 



1 



(ad)r^ sin^ edof'{R), 



-dlf'{R)-3-dof'{R), 



(C.13) 



en donde usamos que el escalar R, NO depende de las coordenadas espaciales, de modo que cualquier 
derivada con respecto a las coordenadas es nula. Asi tendremos que 

[VoVo + □]/'(/?) = -3-ao/(i?), 



[VoVo + n]f\R) = -3-dRf'iR)doR, 



[VoVo + □]/'(/?) = -3-r(i?)i?, 

a 



con R = doR y f"{R) = dfif'{R), con esto la ecuacion de campo nos queda 



d fjR) h f"{R)R ^ Kp 
a^2f'iR) ^ a f'{R) f'{Ry 



sumando y restando el termino 



tendremos 



R fa 



a\ k 

+ 

a 



\ ^ f{R) ^ ^hf"{R)R Kp 



obtenemos finalmente que 



if2 + 



a2 3f'{R) 



a? 2f'{R) a f'{R) f'{Ry 
{RfiR) - f(R)) 



Kp + 



3Hf"{R)R 



(C.14) 

(C.15) 

(C.16) 
(C.17) 

(C.18) 



en donde usamos que -ff = ^ y la expresion para el escalar de curvatura R. Para la componente espacial 
ij tendremos la ecuacion de campo 



f'{R)Rij - \f{R)g^J - [V,V, - gjjW'iR) = KTij, 



(C.19) 
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{ad + 2d2 + 2k)gijf'{R) - -f{R)g^j - [ViVj - g^Ja]f{R) = kq^-TI 
de la expresion (C.7) el termino VjVj es 



(C.20) 



(C.21) 



donde de nuevo usamos que las derivadas espaciales de f'{R) son nulas. Para el termino □/'(i?) ten- 
dremos nuevamente 

(C.22) 



Asf tendremos que 



nf'{R) = -dif{R)-3-dof{R). 



[ViVj - g^J^f'{R) = -amjdof'iR) + 9ijdlf{R) + 3ff.,-9o/'(i?), 

[ViVj - g^p]r{R) = 2aa~gijdof\R) + a^~gijdif'{R), 

[ViV, - g^J^]f'{R) = 2aa~gijf"{R)R + a%^jU"' {R){R? + f"{R)R), 



(C.23) 



de modo que la ecuacion de campo nos queda 
{ad + 2a^ + 2kYgijf\R)-]^a%jf{R)-2aa~gi,f\R)R-a^~g,j{r = ^a^gijT^, (C.24) 



+ 2 



2k 

a J 



f{R) ^af"{R)R r{R){R? f"{R)R_ np 



2f'iR) a f'{R) 



f'{R) 



f'iR) f'iR)' 



sumando y restando el termino 



+ 



a J 



tendremos finalmente 



2H + + 



f'iR) 



Kp + 2HRf"{R) + (/(-^) RfiR)) ^ ^^//(^) ^ {Rff"'{R) 



(C.25) 



(C.26) 



(C.27) 



donde usamos de nuevo la expresion para H. 
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Apendice 



D 



Contribucion de los operadores T>q,q en 
la EDG 



Vamos a considerar la contribucion 

{S^^Vsfs - 62V^p + g^sV^^ - gp-,Vs'')f'{R)V^V'rj\ (D.l) 

usando para esto lo resultados para los operadores usando como espacio-tiempo el de Robertson- 
Walker 

Df{R) = -dlf'iR) - 3Hdof'{R), 

= -f"{R)R - f"\R)R? - mf"{R)R, (D.2) 

Poo = -3Hdof'{R), 

= -3Hf"{R)R, (D.3) 

= 2Hgijdof'{R) + g:3dlf{R), 

= 2Hgijf"iR)R + gij{r\R)R + r{R)R^). (D.4) 

De estas expresiones es claro que los operadores mixtos Pq? = ^io son nulos, y ademas lo son tambien 
contribuciones de la forma Vq y T>^. Reescribamos la expresion (D.l) en la siguiente forma 

+ gpsV.^)nR)V^V'r,^ - {SfV^p + gp^V,^)nR)V^vS^ ■ (D.5) 

Un metodo simple es fijar a para las componentes temporales y espaciales, y determinar la contribucion 
de cada operador. Consideremos entonces cada uno de los dos casos. 



Contribucion a a = 
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Con a = la expresion (D.5) se reduce a 

{5?^Vsi3 + g^sV^')f'{R)V^V'ij^ - {d'^.V.p + gp^V,'') f {R)vPv'r,\ (D.6) 

puesto que las deltas de Kronecker 5^ y los operadores mixtos se hacen cero en ambos parentesis 
para 7 = (5 = i, con i representando las componentes espaciales, la relacion anterior se reduce a 

[Vsp + gpsVo^)f{R)V^V'rf - {V,p + g^.V^^) nR)V^V'^ri\ (D.7) 

expandiendo expli'citamente las sumas, de nuevo recordando que solo los operadores Pqo Y T^ij son 
diferentes de cero tendremos 

[V,p+gp,V^^)nR)V''v'r^'' - {V,p + gp,V,^) nR)V^V^T,^ = 

(Poo + 5oo2?o°)/'(^)^V%o - (Poo + <?oo2?o°)/'(«)^V°r?° + (A,- + g.,Po°)/'(^)^'V^''^° 

- {V,j + giP^^)nR)V'V%\ (D.8) 

{Vsp+gpsV^^)nR)VPv'rf - {V,^ + gfs^V,')f'{R)V^V%'^ = 

(A, + g,,V,'')f'{R)V'V^r,' - {V,, + g,,V,'') f {R)V'V^v' , (D.9) 

escribiendo la forma expHcita de los operadores (D.3) y (D.4), con Vq = g^^V^Q = —Poo tendremos 

{Vsf,+gfssK)f'{R)V''v'7i' - {V,^ + g^^Vo'>)f'{R)V^V%^ = 

{2Hgijf"{R)R + gij{f"{R)R + f"'iR)R^) + +3gi,Hf"iR)R)f'iR)V'V^v" 

- {2Hg,jf"iR)R + g^,{f"{R)R + r{R)R^) + 3gi,Hf"{R)R)f'{R)V'V%^, (D.IO) 



iVsp+gi3s'Do'')nR)V^V'v' - {V,p + g^^V,^) nR)V^V^r,^ = 

g,, {bHf'{R)R + (/" (i?)ij + r{R)R^)+)nR)V'V^r^'' 

-g,,{5Hf"{R)R+{f"{R)R + f"'{R)R'))f'{R)V'V^rf. (D.ll) 

Ahora, escribimos r]aV°' = 

'?a^" = <?oo??V°+5i,??V^=0, =^ 5oor?V° = (D.12) 

de modo que 

{Vsp+gpsV,^)nR)V^V'rf - {V,p + gp,V^^)f{R)V^V^r,^ = 

g,, {5Hf"{R)R + {f"{R)R + f"'{R)R')+)f'{R)V^V^v'' 

+ goo{5Hf"{R)R + {f"{R)R + f"'{R)R'))f'{R)V'V%', (D.13) 
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{5Hf"{R)R + {nR)R + f"{R)R^)+)f{R) [g^.V'V^ + 5ooV^V°r?>°. (D.14) 
Nuevamente la expresion VaV" = e puede expandirse como 

VaV' = 500 + gijV'V^ = e, (D.15) 
de modo que el termino entre parentesis es simplemente e, y asi la contribucion para a = es finalmente 

e{5Hf"{R)R + if"iR)R + r{R)R^)+)fiR)r,'. (D.16) 



{6^^Vs^ + gpsT)^')f'iR)V^V%'' - {6's^,p + g^,V,')f'{R)V''v'rj' 



• Contribucion a a = i 

Con a = i la expresion (D.5) se reduce a 

{Sl^Vsp + g^sV;)nR)V^V'r^'^ - [SlV^p + gp^Vi)f'{R)V^V'r^\ (D.17) 

de nuevo, las deltas de Kronecker 61^ y los operadores mixtos se hacen cero en ambos parentesis 
para 7 = 5 = 0, de modo que 

{Vsp + gps'D;)f'{R)V^V'v' - {V,p + gp.Vl) f {R)VPVW . (D.18) 

aqui es el operador para una sola componente, no debe confundirse con una suma. Expandiendo 
expKcitamente las sumas tendremos 

{Vsp+g(S5V;)nR)V^V'r]^ - {V^p+gp^V:)f{R)VPvW = 

(Poo + 5ooP/)/'(i?)^V^* - (Poo + gooVi)f{R)V^V'if + {V,u + gjkV;)f'{R)V^vW 

- {Vjk + gjkV;)nR)V^V'v\ (D.19) 

{Vsf,+gf,sV;)nR)V^V'rj^ - [V^p + gp^V^) nR)VPvW = 

(Poo + gooVi)f{R)V^V^rf - (Pqo + gmV^) f {R)V^VW (D.20) 

escribiendo la forma de los operadores (D.3) y (D.4), con T>^ = g^^'Dij tendremos 

{Vsp+gp5V,')f'iR)V^V'r]' - {V.,p+gp^Vi)r{R)V^VW = 

{-mf"{R)R - g'^ [2Hg,,f"{R)R + g^,{f"{R)R + f" {R)R')]) f {R)V'^V\' 

- {3HnR)R - 9^^ [2Hg,,f"{R)R + g^,{nR)R + f" iR)R')]) f (RW^V'r,'' , (D.21) 

{Vsf,+gfSsV;)f'iR)V^V'r]' - {V^p+gp^Vi)r{R)V^VW = 

- {5Hf"{R)R + {f"{R)R + f"{R)R^)] ) f'{R)V'vW 

+ {pHf''{R)R + {f''{R)R + f''{R)R?)\)f{R)V'^V'rf, (D.22) 



91 



{Vsp + gfi5Vi)f{R)V^V'r,' - [V^^ + gp^V^) f {R)V^VW = 

goo{5Hf"{R)R + {f"{R)R + f"'{R)R')])f{R)V''V%' 

-goo{5Hf"{R)R+{f"{R)R + f"\R)R')])f'{R)V''V'v", (D.23) 

en donde consideramos que el producto g'^^gij = 1, pues es para una sola componente y usamos que 
500 = —1- Usando de nuevo que goov'^V^ = —gijrfV^ tendremos 

[Vsp + gpsVi)f'{R)V^V'ri' - [V^p + gp^Vi)f'{R)V^VW = 

goo{5Hf"{R)R + {f"{R)R + f'"{R)R')])f'{R)V''V%' 

+ gij{5Hf"{R)R + {f"{R)R + f"{R)R?)])f{R)V^V'r^\ (D.24) 

{Vsp + gpsVi)f\R)V^V'ii^ - {V,p+gp,Vi)nR)V^'VW = 

{5Hf"{R)R + {f"{R)R + r{R)R^)])f'iR) [gooV'V^ + g^JV^Vy'. (D.25) 

Con la expansion para VaV°' = e obtenemos finalmente que la contribucion para una componente r/* es 

{Vsp + gps'Di')f{R)V<'v'r,' - {V^^ + gp^V^')r{R)V^VW = 

e{5Hf{R)R+{f{R)R + r{R)R^)])r{R)r,\ (D.26) 

Asi tenemos entonces que la contribucion es la misma considerando los casos por separado para a = 
y a = i, y que ademas, como era de esperarse, no existen terminos cruzados (con cambio en la direccion 
del vector ry"). Por lo tanto la contribucion total de los operadores es 

{S^Vsfi - S^V^p + gps-D^" - g/s^Vs^)f'{R)V^V'v'' 

= e{5Hf{R)R+{f"{R)R + f"{R)R')])f{R)n-. (D.27) 
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